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ABSTRACT 

This  thesis  examines  the  design  of  the  guidance,  steering,  control,  and 
attitude  measurement  systems  for  the  third  stage  of  a space  vehicle  in 
exoatmospheric  powered  flight.  The  guidance  methods  used  in  this  study 
include  position  offset  guidance  and  energy  management  (fuel  depletion) 
guidance,  davel-epad  hy  The  Chiirlee-Ctarlc  Draper  Laboratory^  A general- 
ized flight  control  system  design  method  is  developed  which  allows  the 
autopilot  filter  gains  and  coefficients  to  be  readily  determined  for 
any  desired  crossover  frequency.  Several  steering  modes  for  velocity- 
to-be-gained  and  fuel  depletion  guidance  are  then  discussed.  The  effect 
of  terminal  steering  on  the  end-of-burn  stability  and  velocity- to-be- 
gained  residuals  is  discussed  for  abrupt  thrust  termination  and  thrust 
tailoff  conditions.  An  empirical  approach  to  the  determination  of 
sensitivities  to  thrust  direction  errors  for  fuel  depletion  and  simul- 
taneous control  of  reentry  angle  and  time-on-target  guidance  is  discus- 
sed and  compared  to  analytical  predictions.  Finally,  an  attitude  meas- 
urement study  is  presented  which  compares  several  methods  of  processing 
the  IMU  attitude  measurements  with  regard  to  reducing  computational 
requirements  and  minimizing  worst-case  error. 
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In  the  design  of  space  vehicle  flight  control  systems  for  highly 
constrained  guidance  methods,  the  system  elements  must  be  carefully 
analyzed  in  terms  of  the  required  vehicle  response  to  guidance  commands. 
These  guidance  commands,  particularly  in  the  case  of  highly  constrained 
methods,  may  require  at  the  same  time  a high  degree  of  steering  accuracy 
and  a well-behaved  response  to  large  transients. 

In  the  case  of  fuel  depletion  steering,  for  example,  if  we  assume 
that  the  vehicle  is  initially  thrusting  along  the  velocity-to-be-gained 
vector,  there  will  be  a large  instantaneous  attitude  error  (50  degrees 
might  be  a typical  value)  as  the  vehicle  maneuvers  to  "waste"  excess 
fuel.  This  large  input  command  can  result  in  extremely  high  turning 
rates,  which,  even  neglecting  the  structural  aspects,  can  severely  tax 
the  capability  of  the  flight  control  system.  This  is  a particularly 
severe  problem  for  flight  control  systems  which  are  designed  to  precisely 
null  the  velocity-to-be-gained. 

The  transient  may  also  cause  saturation  of  the  engine  actuator, 
forcing  the  actuator  into  position  and  rate  limits.  Also,  since  the 
attitude  data  is  sampled  at  discrete  times  and  extrapolated  over  an 
interval,  the  high  turning  rates  can  cause  significant  attitude  (and 
hence  derived  rate)  errors.  The  high  rotation  rates  may  also  cause 
significant  errors  in  the  sensed  acceleration  due  to  the  location  of 
the  IMU  away  from  the  vehicle  center  of  gravity. 

In  spite  of  these  problems,  the  overall  goal  of  the  flight  control 
system  is  to  cause  the  vehicle  to  precisely  achieve  the  velocity  re- 
quired to  coast  to  the  impact  point  at  thrust  termination.  This  re- 
quirement drives  the  design  to  be  highly  sensitive  to  small  steering 
errors  which  is  at  odds  with  the  requirement  to  be  well-behaved  in  the 
presence  of  large  transient  signals. 


In  this  thesis,  these  problems  are  examined  in  terms  of  their 
impact  on  the  design  of  a flight  control  system  for  a Stage  III 
vehicle.  The  system  design  is  evolved  to  minimize  residual  velocity- 
to-be-gained  errors  and  maintain  vehicle  stability  despite  the  highly 
constrained  nature  of  the  guidance  techniques . 
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SECTION  1 


SYSTEM  DESCRIPTION 


1.1  The  Vehicle 

The  vehicle  used  in  this  study  is  the  third  stage  booster  and 
payload  of  a ballistic  missile.  The  missile  is  assumed  to  be  a rigid 
body,  symmetric  about  the  longitudinal  (X)  axis.  The  vehicle  is  powered 
by  a solid  fuel  rocket  engine  with  constant  thrust  (unless  otherwise 
specified)  and  the  vehicle  mass,  moments  of  inertia,  and  center  of  gravity 
are  assumed  to  vary  linearly  with  time.  The  vehicle  is  assumed  to  be 
operating  outside  the  atmosphere.  The  vehicle  parameters  are  illustrated 
in  Figure  1 and  listed  in  Table  1. 
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BODY  COORDINATE  FRAME 
AND  VEHICLE  PARAMETERS 
AND  SIGN  CONVENTION 


Figure  1.  The  vehicle  model 
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Control  of  the  vehicle  is  accomplished  by  deflecting  the  rocket 
nozzle  in  pitch  and  yaw.  The  pitch  and  yaw  autopilot  channels  are 
identical.  Motion  about  the  roll  axis  is  assumed  constrained  to  zero. 

It  is  further  assumed  that  there  are  no  inertial  coupling  effects. 

The  derivations  of  the  vehicle,  actuator,  and  steering  loop 
transfer  functions  are  given  in  Appendix  A.  Although  the  digital  auto- 
pilot operates  by  sampling  inputs  and  performing  calculations  at  dis- 
crete intervals  of  time,  we  will  assume  the  autopilot-actuator-vehicle 
combination  to  be  approximately  continuous,  since  the  sampling  frequency 
of  the  autopilot  is  high  (33.33  cycles/sec)  relative  to  the  expected 
crossover  frequency  (<8  rad/sec) . 

The  open-loop  transfer  function  for  the  actuator  and  vehicle  from 
Appendix  A is  given  below.  "Tail-Wags-Dog"  and  IMU  location  effects 
have  been  analyzed  and  will  be  neglected.  A second-order  actuator  model 
will  be  used.  System  lags,  such  as  sampling  delays  (T/2)  and  computation- 
al lags,  are  combined  into  the  delay  term 


Vehicle  Actuator  Lags 


with 

Kv  = Vehicle  DC  gain  = T£Cg/I 

K = Actuator 
s n 

Tq  = Z delay  terms  (sampling  and  computational) 

The  frequency  response  of  the  uncompensated  vehicle  and  actuator  is 
given  in  Section  2.5,  along  with  the  response  curves  for  several  auto- 
pilot filter  configurations  and  steering  loop  designs. 
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rate,  in  particular,  the  required  velocity  and  the  reentry  angle. 

The  required  velocity,  vreg»  is  the  velocity  needed  to  coast  from  the 
current  position  to  the  desired  impact  point.  The  velocity-to-be-gained, 
Vg,  is  defined  as  the  vector  difference  between  the  required  velocity 
and  the  current  vehicle  velocity 


V 


g 


(1) 


The  reentry  angle  is  the  predicted  flight  path  angle  at  300,000  feet  and 
may  be  controlled  by  selecting  the  appropriate  guidance  mode  (see 
Section  3) . 

The  guidance  parameters  are  computed  by  specifying  the  present 
and  final  position  vectors,  r,  and  ?2,  and  the  desired  time  of  flight 
and  solving  the  "Lambert  Problem. ” The  inertial  reference  frame  and 
the  nominal  trajectory  used  in  this  study  are  illustrated  in  Figure  3. 


Figure  3.  Inertial  coordinate  frame  and  nominal  trajectory. 
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The  flight  control  system  consists  of  the  digital  autopilot, 
steering  loop  and  engine  actuator.  Steering  is  accomplished  at  the 
steering  cycle  rate  by  computing  a desired  thrust  vector  in  platform 
coordinates,  USF,  and  by  taking  the  vector  cross-product  between  the 
incremental  sensed  velocity  vector,  AV,  and  USF.  This  steering  method, 
popularly  known  as  cross-product  steering,  results  in  an  inertially- 
ref erred  steering  vector,  ALPHAS TEER,  where 

ALPHASTEER  - AV  x USF  (2) 

For  simple  cross-product  steering  (no  integral  term  in  the  feedback) 
the  rate  command  vector  in  the  body  axis,  uc,  is  obtained  by  the 
expression 

i oc  = P*B^Kgteer  ALP  HAS  TEE  rJ  (3) 

where  P^  is  the  platform-to-body  transformation  matrix  and  Kgteer  is 
the  steering  gain. 

When  integral  feedback  is  added,  the  body  rate  command  vector  is 
obtained  by  the  operation 

uc  ■ P*B[Ksteer  ALPHASTEER  + Kint  ALPHAINt]  (4) 

where  Kinfc  is  the  integral  steering  gain  and  ALPHAINT  is  computed  each 
steering  cycle  from 

ALPHAINT (nTg)  = ALPHAINT (n-1) Tg  + ALPHASTEER  Tg  (5) 

where  Tg  is  the  sampling  time  of  the  steering  loop. 

The  incremental  attitude  command  vector  to  the  autopilot, 
THCMDINC,  is  computed  by 

THCMDINC  - wcTDAp  (6) 

where  TpAp  is  the  sampling  time  of  the  digital  autopilot. 
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The  nominal  cycle  times  for  the  guidance,  steering,  and  autopilot 
times  are  given  below  in  Table  2. 

Table  2.  Nominal  cycle  times. 


System 

Element 

Cycle  Time 

No.  of  DAP 

Guidance 

<v 

1.35  sec 

45 

Steering 

<TS) 

0.45  sec 

15 

Autopilot 

(tdapj 

0,03  sec 

1 

The  autopilot  computes  an  error  signal  by  comparing  measured  attitude 
increments  with  the  commanded  attitude  increments.  The  compensation 
filter  provides  gain  and  phase  stabilization  to  this  conditionally 
stable  system  by  modifying  the  open-loop  gain  and  adding  phase  lead  to 
meet  the  design  requirements  given  in  Section  2.1.  The  transfer  function 
of  the  compensator  is  given  below 


w + w," 

D(w)  = Kw  w Domain 

2—1 

With  the  substitution  w = — —■ r,  the  function  becomes 

z + 1 

1 — A z"1  1 

2 z Domain 

D(z;  = K zr 

z L 1 - B.z  L J 


where 


- r ^ 1 + w2i  _ 

w2  1 + w3  KW 


1 - v?2 
1 + w2; 


1 - w. 


The  autopilot  DC  gain  is  continually  modified  as  a function  of  the 
magnitude  of  the  thrust  acceleration  in  order  to  maintain  the  open 
loop  gain  at  a constant  value  as  described  in  Section  2.6. 
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, The  Inertial  Measurement  Unit  (IMU)  used  in  this  study,  is  the 

\ Advanced  Inertial  Reference  System  (AIRS)  described  in  detail  in 

Section  4.  This  system  provides  measured  body-angle  increments  as  a 
feedback  signal  to  the  autopilot  loop  and  accelerometer  (assumed  ideal) 
sensed  velocity  increments  as  a feedback  signal  to  the  guidance  and 
steering  systems.  The  attitude  feedback  approach  is  to  initialize  the 
attitude  errors  at  the  beginning  of  the  burn  and  then  to  update  the 
attitude  errors  approximately  by  adding  to  the  attitude  error  for  each 
body  axis  the  difference  between  the  commanded  body-angle  increment 
and  the  AIRS  measured  body-angle  increment  for  that  axis  every  auto- 
pilot cycle.  The  advantages  of  this  approach  are  (1)  it  requires  only 
one  time-consuming  computation  of  the  platf orm-to-body  transformation 
matrix  from  measured  band  angles  for  error  initialization  at  the  be- 
ginning of  the  burn,  and  (2)  it  also  minimizes  the  computation  time  by 
| precisely  undating  the  band-angle-increment-to  body-angle-increment 

ij  matrix  only  every  steering  cycle. 

i 

tj  After  initialization,  computation  time  is  minimized  by  using 

| feedback  body-angle  increments  determined  by  multiplying  the  measured 

| band  angle  increments  over  autopilot  cycle  by  a band-angle-to-body-angle 

j transformation  matrix.  This  incremental  transformation  matrix  requires 

» less  computation  time  than  computing  the  platform-to-body  transformation 

I matrix,  and  does  not  have  to  be  updated  precisely  every  autopilot  cycle. 


r t. 
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SECTION  2 


STEERING  AND  AUTOPILOT  ANALYSIS 


2.1  Steering  and  Autopilot  Design  Parameters 


The  overall  function  of  the  steering  and  autopilot  combination 
is  to  minimize  steering  errors  without  exceeding  actuator  position  or 
rate  limits.  The  steering  errors  are  generally  reflected  by  velocity- 
to-be-gained  errors  at  thrust  termination.  Another  requirement  is  that 
the  vehicle  response  be  stable  throughout  the  burn  and  particularly 
near  the  end,  when  the  steering  methods  used  in  this  study  are  inherently 
unstable.  As  we  shall  see,  the  steering  mode  must  be  switched  to  a 
"terminal"  mode  near  the  end  of  the  burn,  and  it  may  also  be  advantageous 
to  vary  the  autopilot  parameters  at  this  time. 

In  terms  of  specific  autopilot  design  goals,  we  shall  specify 
the  phase  margin,  gain  margin,  and  peak  closed-loop  gain.  In  addition 
we  shall  also  specify  residual  Vg  "goals".  These  design  paramters  are 
given  in  Table  3 below. 


Table  3.  Design  parameters. 


Parameter 

Open-Loop  Gain  Margin 
Open-Loop  Phase  Margin 
Peak  Closed-Loop  Gain 
Residual  Velocity  to  be  Gained 


Steering  Modes 


Design  Value 


30  degrees 


<10  fps 


The  vehicle  may  be  operated  in  either  of  two  distinct  steering 
modes;  Vg  (velocity-to-be-gained)  steering,  or  fuel  depletion  steering. 
In  Vg  steering,  the  vehicle  thrust  vector  is  pointed  in  the  direction 
of  the  velocity-to-be-gained  vector,  which  is  computed  by  the  guidance 
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computer  as  the  difference  between  the  required  velocity,  vreq'  an<*  the 
present  velocity,  V.  In  vector  notation,  this  is  again 


\7  = V 

g req 


(1) 

(Repeated) 


The  required  velocity  at  a given  point  on  the  trajectory  is  calculated 
by  solving  the  Lambert  problem,  given  the  present  and  desired  final 
radius  vectors,  and  the  desired  time  of  flight.  The  vehicle  radius 
vector  is  offset  to  account  for  the  nonimpulsive  nature  of  the  thrust. 
This  guidance  scheme  is  discussed  in  detail  in  Section  3.  steering 

assumes  that  the  vehicle  has  thrust  cut-off  capability  or  that  the 
thrust  may  be  modulated  so  that  the  AV  capability  of  the  vehicle  is 
exactly  equal  to  the  magnitude  of  the  velocity-to-be-gained.  In  V 
steering,  then,  the  vehicle  is  steered  such  that  the  thrust  vector  is 
pointed  along  the  velocity-to-be-gained  vector  and  the  engine  is  shut 
off  when  the  magnitude  of  the  velocity-to-be-gained  goes  to  zero.  Since 
all  three  components  of  the  Vg  vector  do  not  normally  go  to  zero  at 
exactly  the  same  time,  the  engine  is  usually  shut  down  when  the 
component  of  along  the  vehicle  x-axis  goes  to  zero. 

If  the  vehicle  has  a AV  capability  which  is  greater  than  the 
magnitude  of  the  velocity-to-be-gained  for  a given  burn,  and  the  engine 
cannot  be  shut  off,  then  fuel  depletion  steering  must  be  used.  In 
this  steering  mode,  the  thrust  vector  is  rotated  at  an  angle  to  the  Vg 
vector  to  "waste"  the  excess  fuel.  As  the  magnitude  of  the  AV  capability 
is  reduced,  the  angle  between  the  thrust  vector  and  the  velocity-to-be- 
gained  vector  goes  to  zero  in  such  a way  that  when  the  fuel  depletion 
angle  is  zero,  the  remaining  AV  capability  and  the  velocity-to-be-gained 
are  also  zero.  The  geometrical  relationships  of  fuel  depletion  steering 
are  shown  in  Figure  4.  After  the  vehicle  rotates  to  the  initial  fuel 
depletion  angle,  the  rate  of  change  of  the  fuel  depletion  angle  is 
almost  constant.  As  shown  in  Figure  4,  the  vehicle  thrust  vector 
rotates  through  an  angle  26  during  the  burn,  so  this  rate  is  approx- 
imately 26/Tburn. 

When  fuel  depletion  steering  is  initiated,  a plane  containing  the 
radius  vector  r,  the  Vg  vector,  and  the  normal  to  the  Vg  vector,  the 
unit  vector  B,  is  defined  in  velocity  space.  The  fuel  depletion  angle 
6 is  computed  from  Equation  (7)  which  is  derived  in  Appendix  4 


sin  6 


26 
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Figure  4.  Geometry  of  the  fuel  depletion  arc.  \ 

This  relationship  is  obtained  from  the  geometry  of  the  arc  of  a circle 
where  20  is  the  interior  angle  subtended  by  a chord  line,  represented 
by  Vg,  and  the  arc  length  is  equal  to  the  magnitude  of  the  AV  capability. 

A } 

The  desired  thrust  vector,  USF,  is  computed  from  0,  Vg,  and  B 
using  the  equation  < 

USF  = cos  0 Unit  (V  ) + sin  0 (B)  (8) 

9 


< 


The  USF  vector  is  always  tangent  to  the  circular  arc  and  forms  the 
angle  0 with  the  Vg  vector. 

The  computation  of  the  fuel  depletion  angle  0 may  be  simplified 

by  expanding  (sin  0)/0  in  a power  series  and  neglecting  terms  of  order 

4 

0 . This  gives  the  simpler  expression  given  in  Equation  (9) 

8 = '6(1  - Vg/Vcap)  (9) 

A comparison  of  the  resulting  0 using  this  method  and  the  exact 
method  and  the  sensitivity  of  the  fuel  depletion  angle  to  small  Vcap 
changes  are  given  in  Table  4.  Note  that  the  angle  computed  by  Equation 
(9)  is  always  less  than  the  exact  method. 
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Table  4.  Sensitivity  of  the  fuel  depletion  angle  6 

to  a small  change  in  the  ratio  of  V /V 

g cap 


V /V 
g cap 

0 (deg) 
exact 

A0  (deg) 

% Change 

8 (deg) 
simplified 

+2% 

, 0.816 

61.99 

-2.82 

-4.35 

Nominal  , 

t 0.800 

64.81 

- 

- 

62.76 

-2% 

( 0.784 

67.54 

2.73 

4.21 

- 

+2% 

0.918 

40.70 

-4.37 

-9.70 

_ 

Nominal  | 

1 0.900 

45.07 

- 

- 

44.38 

-2%  1 

! 0.882 

49.10 

4.03 

8.94 

- 

+2%  - 

0.969 

24.83 

-6.79 

-21.47 

_ 

Nominal  j 

1 0.950 

31.62 

- 

- 

31.38 

-2%  1 

! 0.9310 

37.26 

5.64 

17.84 

- 

+2% 

. 0.997 

2.81 

-17.10 

-85.89 

Nominal 

i 0.980 

k 

19.91 

- 

- 

19.85 

-2%  1 

l 0.960 

28.10 

8.19 

41.14 

_ 

2.3  Steering  Stability  Analysis 

The  stability  characteristics  of  the  steering  modes  presented  in 
this  study  are  conveniently  analysed  in  terms  of  the  frequency  response 
characteristics  of  the  pitch  (or  yaw)  autopilot  channel  and  steering 
loops.  As  previously  stated,  we  assume  that  the  roll  orientation  is 
constant,  and  translational  motion  occurs  in  the  pitch  (or  yaw)  plane. 

The  single  plane  analytical  model  to  be  used  in  this  analysis  is 
presented  in  Figure  5.  This  model  represents  the  steering  effects  in 
terms  of  two  sampled-data  transfer  functions  F-^*  (s)  and  FjMs)  which 
are  based  on;  the  steering  sample  period  Tg.  The  key  definitions, 
assumptions,  and  approximations  of  this  model  are: 

(1)  The  autopilot  and  steering  systems  are  approximated  by 

linear  perturbation  equations  where  the  vehicle  rotation, 

9,  the  engine  deflection  6,  and  the  inertial  rotation  of 
the  thrust  vector,  Ay,  have  zero  steady-state  values.  The 
steady-state  position  of  the  vehicle  center  of  mass  is 
assumed  to  be  on  the  vehicle  x-axis,  and  the  steady-state 
directions  of  the  commanded  and  actual  thrust  vectors  are 
assumed  to  be  along  the  steady-state  inertial  orientation 
of  the  vehicle  x-axis. 
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AUTOPI  LOT  LOOP  APPROXIMATE  CONTINUOUS  MODEL 


CROSS  PROOUCT  STEERING  BASED  ON  AV  x USF 
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T - 1 Vg  STEERING 
7*4  FUEL  DEPL  STEERING 


INTEGRAL  STEERING 
STEER02 
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*USF*‘>  9°  I 

Figure  5.  Autopilot  and  steering  loop  analytical  model. 

(2)  Since  the  autopilot  sampling  period  TDAp  is  much  shorter 
than  the  steering  period  Tg,  it  is  possible  to  treat  the 
autopilot  as  continuous  and  to  combine  the  autopilot  func- 
tion 0 (s)/ec(s)  with  other  factors  to  obtain 


Fl(s) J 


1 / 6 (s)  \/  X (s)  \ 

3(ec(s)^e(s)j 


(3)  It  can  be  shown  that  the  sampled  function  F^*(s)  which  is 
expressed  as 


oo 

F1MS)  = T'  ^ F1.(S  + 

8 n»-°° 

can  be  approximated  by  the  n=0  term 


Fl*(n)  = ah  F,(S) 
s 


■r  % • i. — t rr 


and  so 


F1*(s) 


(mtm 


This  approximation  can  be  justified  on  the  basis  that  the 
1/s3  term  in  F^(s)  causes  the  |n|>0  terms  in  F1*(ju)  to 
be  small  compared  to  the  n=0  term  in  the  frequency  range 
less  than  half  the  steering  sampling  frequency. 

Starting  with  the  approximation  that  F^Ms)  = 1/Tg  F^(s) 
and  noting  that  the  resulting  steering  loop  transfer  function 
is  given  by 


Fx*(s)  F2*(s)  = 


F2* (3) 


it  can  be  reasoned  that  the  effect  of  this  transfer  function 
can  be  equivalently  represented  by  treating  the  steering  as 
merely  an  addition  to  the  unity-gain  attitude  feedback  of 
the  autopilot.  This  yields  a modified  "total"  open-loop 
function  of  the  autopilot  which  is  expressed  as 


,n>  _ [ 6(8)1  . / 1 \i  1 \/X(s)V 

(8)  " |^j  1 


f2*<s) 


(5)  The  sampled-data  transfer  function  F2*(s)  is  expressed,  as 
shown  in  Figure  5,  in  terms  of  the  sampled  integral  Xy  of 
the  thrust  angle  X and  in  terms  of  the  sensed-velocity  angle 
X^v  and  commanded  thrust  angle  ^USF- 

(6)  The  sensitivity  of  the  angle  X„__,  of  the  commanded  thrust 
direction  to  the  sampled  angle  Xv  of  the  integral  of  the 
thrust  angle  X is  approximately  expressed  as  the  ratio 
y/TgQ,  where  y=l  for  Vg  steering  and  y=4  for  fuel  depletion 
steering,  and  where  Tgo  is  the  time-to-go.  The  derivation 
of  these  sensitivities  is  presented  in  Appendix  A4  and 
Section  3.3. 

(7)  An  approximate  expression  for  X(s)/0(s)  is  given  in  Eq.  (A-12) 
of  Appendix  A3.  It  can  be  shown  that  for  the  vehicle  para- 
meters used  in  this  study  that  X(s)/0(s)  = 1,  and  this  assump- 
tion will  be  made  for  the  purposes  of  the  following  analysis. 


He  now  consider  two  steering  models,  STEER61  (shown  in  Figure  6) 
in  which  the  gain  of  the  integral  path  of  the  steering  law  is  set  to 
zero  and  STEERS 2 (shown  in  Figure  7)  which  includes  the  integral  path. 
In  case  STEERS  1,  the  function  F2*  (s)  is  given  by 


F2*(s) 


_Ksteer^1 


- V\  1 - e*  V 


Assume  for  this  case  that  T is  sufficiently  large  such  that  y/T  can 

go  go 

be  neglected 


V<s)  5 -Ksteer 


/ 

Now,  I 1 - e ] can  be  transformed  first  to  the  z-domain-and  then  to 
the  w-domain  by  the  substitutions 


-T  s -1 

e = z 


_ 1 - w 

1 + w 


Therefore 


F2*(s)  = 


4 /_w\2 
Tg  (l  + 


w is  related  to  real  frequency,  w,  by  the  relationship 


w = j tan 


fr) 


Assume  that  the  steering  loop  crossover  frequency  is  much  smaller  than 
the  half-steering  loop  sampling  frequency.  Then  for  the  region  of 
frequencies  of  interest 


mt  . 25  a 

w - — 


rhen , F2*  (s)  can  be  approximated  as 
F2*(s)  = -K 


steer  s (1  + Tgs) 


The  modified  total  open-loop  function  then  becomes 

1 


TOT 


'd(s)' 

EXIT 


1 4.  ^ (S) 

9 (s)  *steer  s (1  + T s) 


If  we  make  the  further  approximations  that  Tgs  <<  1 and 


then 


Ms) 
0 (s) 


G * r e (s)  1 

, , Ksteer 

6 (s)  s + Ksteer 

^TOT  ITiTl 

1 + s 

Eli!  s 

(10) 


This  simplified  steering  loop  model  is  used  in  Section  2.4  to  select 
the  autopilot  parameters  in  the  normalized  design  process.  The  selec- 
tion of  the  value  of  Ksteer  is  dependent  on  the  design  of  the  autopilot 

and  its  desired  crossover  frequency,  since  we  desire  the  frequency  of 
the  steering  loop  to  be  significantly  lower  than  the  autopilot  crossover 
frequency. 


(«) 


r - 1 FOR  VG  STEERING 
y - 4 FOR  FUEL  DEPLETION  STEERING 

Figure  6.  Analytical  steering  model  - case  STEER01 . 
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In  the  second  steering  model,  (STEER02)  to  be  considered,  the 
effect  of  adding  integral  feedback  to  the  proportional  loop  will  be 
considered.  The  integral  feedback  is  added  to  reduce  the  steady-state 
errors  due  to  the  approximately  constant  rate  command  during  fuel 
depletion  steering.  It  will  be  shown  how  the  integral  gain,  Kint»  can 
be  selected  based  on  Ksteer-  The  STEER02  model  is  described  in  Figure  7 
For  this  case  the  function  F2*{s)  is  given  by 


F2Ms) 


* Ksteer  + 


Kint  Ts 


1 - e 


1 - e 


s°  1 - e 


7 - 1 FOR  VG  STEERING 

7-4  FOR  FUEL  DEPLETION  STEERING 

Figure  7.  Analytical  steering  model  - Case  STEER02. 


If  the  same  assumptions  are  made  as  in  the  case  of  STEER01,  then 

/ -T  s\ 

^1  - e J is  again  replaced  by  sTg  and  y/^g0  is  assumed  negligible. 

THE  STEER02  steering  model  of  Figure  7 reduces  to  that  shown  in  Figure  8 
and  F2* (s)  is  given  by 


F2*(s)  = -sTg 


(Ksteer  3 + Kint^ 


Figure  8.  Simplified  analytical  model  STEER02. 

For  the  assumption  that  A(s)/6(s)  = 1,  the  modified  total  open-loop 
function  (for  which  the  unity-gain  attitude  feedback  is  added  to  the 
steering  loop  feedback)  then  becomes 


gtot(s)  = [Ufr] 

. Ksteer  s + Kint 

s2 

L J 

■ [«»] 


s + Ksteer  s + Kint 


(11) 


If  the  second  order  numerator  is  expressed  in  the  "standard"  second 
order  form 


then 


Therefore 


+ Ksteer  3 + Kint 


2 2 
s + 2£oo  s + w 
n n 


steer 


4^“n  " 4^  Kint 


Select  a damping  ratio  for  the  steering  loop  of  c 
Then 

2 

K . . = K \ 

int  steer 


= 0.5. 
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This  relationship  will  be  used  to  determine  K^nt  after  Ksteer  has  been 


determined  from  the  normalized  compensation  design. 


Finally,  STEER02  can  be  modified  to  include  the  case  near  the  end 


of  the  burn  when  TgQ  becomes  small.  In  this  case,  the  function  FjMs) 
is  approximated  by 


F2*(S)  = -Ts 


[Ksteer  8 + Kint  8 + 

L J g° 


Therefore 


V B2  + l K3fceer  X + V | „ . Y Kint 

Ksteer  8 IT  Kint ) 8 + T 

= lUfjJ  1 * 1—1°.  / £L_ 
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Approximations  for  the  modified  total  open-loop  steering  and 
control  transfer  function  have  been  derived  for  three  different  sets 
of  assumptions  for  the  steering  loop.  These  are  summarized  as  follows: 


Assumptions  for  all  3 cases: 


(1)  For  frequency  range  below  steering  loop  crossover  frequency, 


the  expression  1 1 - e I can  be  approximated  by  T 


(2)  X(s)  - , 

57s) 


Case  I:  steering  Model  STEER01 — No  integral  path  in  steering  law. 


Assumption:  TgQ  large  ^ 


negligible^ 


GTOT 


9 (s)  (s  + KsteerJ 
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Case  II:  Steering  Model  STEER02 — Proportional  and  integral  paths  in 
steering  law. 

Assumption:  T^q  large  negligible  ^ 


,«■>  ■ m[— 


Ksteer  8 + Kint 


(K^nt  is  chosen  to  be  equal  to  Kgteer  corresponding  to  a damping  ratio 
of  0.5) 

Case  III:  Steering  Model  STEER02 — Proportional  and  integral  paths  in 
steering  law. 


Assumption:  T small  I =-*—  term  cannot  be  neglected' 

9 V 9°  I 
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2.4  Normalized  Compensation  Design 

In  this  section,  we  develop  a design  approach  that  is  based  on 
normalizing  the  system  parameters  such  that  the  characteristics  of 
the  autopilot  and  steering  loop  become  independent  of  the  autopilot 
crossover  frequency.  The  result  is  a general  expression  for  the 
simplified  open-loop  transfer  function  that  is  used  to  select  values 
of  normalized  gain  and  system  delays  that  provide  the  desired  response. 
After  these  values  are  selected,  the  system  parameters  are  readily 
determined  for  any  desired  value  of  autopilot  crossover  frequency. 

To  develop  this  expression,  we  write  the  open-loop  transfer 
function  of  the  combined  compensator,  actuator  and  vehicle  derived 
in  Section  1.1  as 


G(s)  = 


<(s/«2)  + l)KcKve 


-Vs 


Equation  (13)  is  obtained  by  treating  the  pole  of  the  rigid  body  compen- 
sator, oij,  and  the  actuator  transfer  function  as  pure  delays  since  they 


1 


occur  at  high  frequencies  relative  to  the  desired  range  of  crossover 
frequencies.  These  lags  are  lumped  into  T ' along  with  the  sampling 
and  computational  delays.  If  we  then  let 


sA,2,  Ko* 


we  can  rewrite  (13)  as 


G'  (s) 


(s'  + 1)  K^e 


We  then  obtain  the  open  and  closed-loop  frequency  response  of  (14)  and 
then  plot  the  locus  of  maximum  values  of  the  closed-loop  gain  H,  M 
against  Kq'  and  <t>0'.  This  data  is  tabulated  in  Table  5 and  plotted  in 
Figure  9 for  autopilot  loop  (without  steering) . Figure  9 indicates 
that,  in  order  to  keep  Mp  less  than  6 dB,  for  = 0.30,  the  allowable 

range  of  Kq'  is  approximately  1.2  £ K0’  < 1.75.  For  4>0'  = 0.25,  KQ ' 


has  an  allowable  range  of  0.8  < Kq' 


If  we  select  a nominal  Kq'  = 1.4, 


for  ^o'  = 0-30'  the  closed-loop  gain  margin  is  0.12  dB  and  for  <j>o ■ = 0.2 
the  closed-loop  gain  margin  is  1.03  dB. 


We  now  turn  to  the  design  of  the  system  including  steering  effects 
For  the  purpose  of  this  design  process,  we  will  use  the  simple  cross- 
product  steering  transfer  function  developed  in  Section  2.3.  The  in- 
corporation of  this  simplified  model  gives  an  open-loop  transfer 
function  shown  in  Equation  (15). 


(s/w2  + DKcKve 
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o 


Again,  we  normalize  this  function  to  the  frequency  of  the  rigid  body 
compensator  zero,  u>2  , and  the  result  is 


r s'  + ml  1 L*  + l\K _'e  4,0  8 1 


rr 


where 


S/(D- 


W1  ‘ = “]/“2 


K ' = 

o 


K K 
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As  in  the  case  without  steering,  this  design  will  result  in  a 0 dB 
crossover  frequency  which  is  proportional  to  w2*  Again,  we  obtain  the 
open  and  closed-loop  frequency  response  of  (16)  and  plot  the  locus  of 
peak  closed-loop  gain,  Mp  versus  Kq'  and  <P0‘ • These  results  are  tab- 
ulated in  Table  6 and  graphed  in  Figure  10,  and  show  that,  with  steering 

included,  the  maximum  value  of  d>  * that  gives  M < 6 dB  is  now  0.25. 

o p — 

For  Kq'  = 1.4,  the  closed-loop  gain  margin  is  0.10  dB  for  $ ' = 0.25 
and  0.91  dB  for  $ ' = 0.20. 

If  we  now  select  a ratio  of  the  steering  frequency  to  the  frequency 
of  the  compensator  zero,  u^/a >2,  of  0.10,  for  example,  we  can  determine 
the  approximate  crossover  frequency  of  the  "real"  transfer  function  as 
a function  of  the  normalizing  frequency  to2- 

Equation  (16)  is  written  as 


G'  (s) 


Ko'  (S'  + iql/u>2)  (s'  + l)e 
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(17) 


With  <**j/o)2  = 0.1,  we  can  safely  neglect  the  gain  contribution  of  the 
steering  loop  at  the  expected  0 dB  crossover  frequency  of  the  autopilot 
loop  and  set  | G * ( jo»)  | = 1. 


K ' (s’  + l)e 

G (s)  = -2 6 

° ..2 


a 


and 


which  is  written  as 


The  positive  solution 
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: | Now  w*  can  be  solved 

| | Table  7. 

for  in  terms 

of  K ' . 
o 

These  results  are 

shown 

F 

? 

Table  7. 

K ' vs 
o 

u'  . 

o 

00 

o 

0 

1.0  1.2 

1.4 

1.6  1.8 

2.0 

i 

; to  * 1.09 

1.27  1.46 

1.64 

1.82  2.01 

2.20 

1 m*4+o>'2  + 


,2 


1 = 0 


1 2 2 2 
■^(0)^)  + (<o'z)  +1  = 0 


for  w'  is 


K «2  /K  '4  ,\1/2 

-v  + (-v  + v) 


1/2 


(19) 


The  crossover  frequency,  uxo*  is  given  by 


10  — Cu  D9 

xo  * 


(20) 


If  we  use  the  previously  selected  value  of  K ' = 1.4,  this  gives  to'  = 1.64, 
and  hence,  is  determined  for  any  desired  crossover  frequency  (approxi- 
mate) . Once  0)2  is  determined,  all  other  system  parameters  are  determined 
from  Equation  (16) . 

The  determination  of  these  specific  parameters  for  autopilot 
crossover  frequencies  of  2,  5,  and  8 rad/sec  is  given  in  Appendix  A5. 

The  results  of  these  calculations  are  summarized  in  Tables  8,  9,  and  10. 
From  Table  8,  we  see  that,  if  we  restrict  the  lead  ratio,  to 

values  less  than  15,  then  a crossover  frequency  of  10  rad/sec  is  not 
attainable. 
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Desired 

Crossover 

Frequency 


Table  8.  Autopilot  parameters. 

s-Plane  w-Plane 


w-Ptane 

Filter 

Gain 


wx0 

H 

B 

B 

W1 

w2 

B 

**K 

*c 

2 

0.12195 

1.2195 

5.9113 

0.001829 

0.01829 

0.08867 

0.10410 

5 

0.30488 

3.0488 

21.6606 

0.004573 

0.04573 

0.3249 

5* 

0.60000 

3.0488 

15.000 

0.09000 

0.04573 

0.22500 

0.65066 

8 

0.48781 

4.87805 

64.8635 

0.007317 

0.073171 

0.97295 

1.66568 

10 

0.60976 

6.09756 

193.548 

0.009146 

0.091463 

2.9032 

am 

*ForcJx0  = 10;  = 31.74  >15  (We  assume  that  this  lead  ratio  is  too  large  to  be 

“2  obtained  with  a simple  compensator.) 

**For  Kv  = 20 
^Optimized  design 


Table  9.  Z-domai’n  filter  coefficients.  j 


DESIRED  wx0 

A 2 

b2 

Kz 

2.0 

0.964072 

0.837104 

0.471993 

5.0 

0.912540 

0.509548 

3.6485 

5.0* 

0.912540 

0.632651 

2.732 

8.0 

0.866636 

0.013710 

12.0475 

^Optimized  design 
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Table  10.  Comparison  of  design  results. 


Desired  u>xq 
(rad/sec) 

Actual  to. 

xo 

(rad/sec) 

Gain  Margin 

(dB) 

Phase  Margin 

(deg) 

Freq.  at  Phase 
Crossover 
(rad/ sec) 

2.0 

1.9 

20.12 

34.6 

9.4 

5.0 

4.9 

12.59 

32.5 

15.6 

5.0* 

4.8 

12.70 

27.0 

13.3 

8.0 

m 

9.2 

30.6 

20.9 

’Optimized  design 

For  the  specific  point  design  to  be  used  in  this  study,  an  autopilot 
crossover  frequency  of  5 rad/sec  will  be  selected  and  evaluated.  If 
this  design  can  be  shown  to  satisfy  the  responsiveness  required  by 
the  steering  techniques  used,  it  will  also  have  the  additional  benefit 
of  reduced  sensitivity  to  noise  from  the  IMU  relative  to  the  8 rad/sec 
autopilot. 

2.5  Frequency  Response 

This  section  contains  the  frequency  response  plots  of  the  uncom- 
pensated system,  (Figures  11-13)  the  compensated  system  resulting  from 
the  normalized  compensator  design,  (Figures  14-22)  and  various  steering 
parameter  effects.  Also,  the  frequency  response  of  the  compensated 
system  for  the  5 rad/sec  autopilot  is  "optimized"  to  obtain  a desired 
closed  loop  gain  over  the  maximum  possible  range  of  frequencies.  Unless 
otherwise  stated  in  these  analyses,  the  nominal  value  of  the  steering 


loop  parameters  Ksteer  (steering  gain),  K^nt  (integral  gain),  TgQ  (time 

to  go  until  the  end  of  the  burn) , and  T (sample  period  of  the  steering 

s 

. loop)  are  used.  These  are  given  in  Table  11. 

Table  11.  Nominal  steering  loop  parameters. 
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The  frequency  response  of  the  steering  models  is  obtained  by  using  the 
exact  steering  representation  from  Figure  5 and  substituting 


~T  s 

1 - e s I = 1 - {cos  wT_  - j sin  wT  ) 

' s a j«»  s s 


If  we  now  say 


C = 1 - cos  uT„ 

U)  5 


s ■ sin  uT„ 

0>  a 


The  total  steering  loop  can  be  represented  as 


Psteer(s> 


* f1  < 


Equation  (21)  is  used  in  determining  the  frequency  responses  given  in 
Figures  23  through  32. 

Figure  23  shows  the  frequency  response  of  the  combined  system 
with  no  steering  and  with  Vg  steering  with  and  without  integral  com- 
pensation. The  other  steering  loop  parameters  are  set  to  their  nominal 
values.  Note  that  in  this  frequency  range  there  is  very  little  dif- 
ference in  the  two  steering  modes  and  both  modes  contribute  phase  lead 
at  the  autopilot  crossover  frequency.  Figures  24  through  29  show  the 
effects  of  changing  the  steering  loop  parameters  Tg0,  Ts)  on 

the  frequency  response  of  the  steering  loop  alone.  Figure  24  shows  that 
increasing  K^nt  gives  increased  gain  at  very  low  frequencies,  but  more 
attenuation  at  about  1 rad/ sec.  Figure  25  shows  that  increasing  t 
gives  more  phase  lag  at  very  low  frequency,  and  slightly  more  phase  lead 
near  the  autopilot  crossover.  The  effects  of  decreasing  the  time-to-go 
are  shown  in  Figures  26  and  27.  The  low  frequency  gain  and  phase  lag 
from  steering  are  increased  significantly  as  TgQ->-0.  Figures  28  and  29 
show  the  effects  of  sampling  period  T for  both  steering  modes.  At 
low  frequency,  the  effect  of  changing  the  sample  period  is  almost 
negligible,  and  at  the  higher  frequencies  where  it  would  become  important, 
the  effect  of  the  other  steering  parameters  causes  the  total  effect  of 
the  steering  to  disappear. 


Figure  24.  Effect  of  gain  vs  frequency 


Figure  26.  Effect  of  T , gain  vs  frequency. 


FREQUENCY  (rad/tec) 


(dB> 


OPEN-LOOP  PHASE  (dig) 

-210  -200  -190  -180  -170  -160  -160  -140  -130  -120 
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OPEN-LOOP  PHASE  (deg) 


This  series  of  plots  show  that  for  = 0.10,  the  steering 

effects  have  little  impact  on  the  autopilot  in  the  vicinity  of  crossover. 
We  would  like  to  use  the  response  characteristics  of  the  steering  to 
modify  the  overall  system  characteristics  to  shape  the  gain-phase  curve 
to  the  M = 6 dB  contour  to  give  an  "optimum"  closed-loop  gain  over  a 
wider  range  of  frequencies.  To  do  this,  we  increase  the  frequency  of 
the  steering  loop  by  changing  the  steering  gain  from  Ksfceer  * 0.3048 
to  Ksteer  = 0.60.  This  corresponds  to  a ratio  of  w^/iOj  =0.20.  We 
also  wish  to  shape  the  gain-phase  curve  at  frequencies  at  and  above  the 
crossover  point.  We  do  this  by  decreasing  the  frequency  of  the  com- 
pensator pole  from  = 21.66  rad/sec  to  u>3  = 15.0  rad/sec.  This  will 
reduce  the  overall  phase  lead  of  the  compensator  and  bend  the  high 
frequency  part  of  the  gain-phase  plot  around  the  6 dB  contour.  The 
effect  of  the  steering  parameters  on  the  steering  loop  response  will 
remain  the  same  except  for  the  frequency  shift,  but  the  steering  loop 
response  now  will  have  more  effect  on  the  total  system  response. 

In  this  optimization  process,  we  use  the  time-to-go  sensitivity 
of  the  fuel  depletion  steering  as  a worst  case,  since  we  have  seen  how 
the  time-to-go  contributes  significant  phase  lag,  and  the  time-to-go 
sensitivity  of  the  fuel  depletion  steering  is  4 times  as  great  as  Vg 
steering.  The  effect  of  steering  on  the  "optimized"  autopilot  is 
shown  in  Figure  30.  Since  the  autopilot  is  optimized  around  the  effect 
of  the  steering  loop,  in  particular  the  phase  lead  contribution,  with- 
out steering  the  gain-phase  curve  crosses  the  6 dB  contour.  The  integral 
term  in  STEER02  gives  more  phase  lag  at  low  frequency,  but  more  phase 
lead  near  crossover.  The  gain-phase  curve  for  steering  mode  STEER01 
crosses  the  6 dB  contour  between  3 and  5 rad/sec,  but  with  a small 
increase  in  open-loop  gain,  this  curve  could  be  adjusted  to  be  tangent 
to  the  6 dB  contour. 

Figure  31  shows  the  effect  of  decreasing  time-to-go  for  STEER01. 
Note  the  dramatic  change  in  the  gain-phase  curve  for  T^o  = 2 sec.  Here, 
the  curve  crosses  the  6 dB  contour  between  1.4  and  .6  rad/sec,  but  this 
can  easily  be  prevented  by  a slight  change  in  open-loop  gain. 

Figure  32  shows  the  effect  of  decreasing  time-to-go  for  STEER02 . 
Here,  the  phase  lag  contributions  from  the  integral  term  and  from  time- 
to-go  cause  a significant  invasion  of  the  6 dB  contour  which  can  only 
be  prevented  by  "freezing"  the  steering  command  at  T values  greater 
than  2 seconds.  In  this  case,  the  dependance  of  the  steering  loop 
on  time-to-go  disappears  entirely. 
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2.6  Gain  Scheduling 

The  analysis  and  parametiic  design  of  the  autopilot  and  steering 
loops  assume  that  the  open  loop  gain,  KCKV  of  the  compensator  and 
vehicle  combination  is  constant.  Since  the  vehicle  gain,  Kv,  is  defined 
as  Kv  - TlCg/1  and  these  parameters  vary  as  the  burn  progresses,  the 
vehicle  gain  increases  from  Kv  = 20  at  ignition  to  * 80  at  burnout 
(assuming  constant  thrust) . When  the  thrust  tailoff  model  is  used, 
the  vehicle  gain  changes  from  a maximum  value  of  Ky  * 61.4  to  K?  =>  0. 

In  order  to  present  the  dynamics  of  the  system  from  drastically  changing 
due  to  these  gain  shifts,  it  is  desirable  to  schedule  the  generalized 
compensator  gain  Kc  as  an  inverse  function  of  the  vehicle  gain  to 
preserve  the  desired  value  of  open-loop  gain.  This  is  accomplished  by 
computing  the  digital  filter  gain  as  follows 

K2  = KV/Kv 

KV  is  given  as  a function  of  the  autopilot  configuration  of  the  initial 
values  of  Kw  (w-plane  equivalent  gain)  and  Kz  (digital  filter  equivalent 
gain)  in  Table  12.  When  the  thrust  tailoff  model  is  used,  the  compen- 
sator gain  is  increased  as  the  vehicle  gain  is  reduced.  In  order  to 
prevent  the  compensator  gain  from  becoming  infinite,  it  is  frozen  when 
it  reaches  its  initial  value.  This  occurs  at  approximately  2 seconds  to 
go. 


Table  12.  Gain  scheduling  autopilot  coefficients. 


Desired  Wxq 

2 

5 

8 


Initial  K 

w 

0.104103 

0.6S0663 

1.665679 


Initial  K 

0.47199 

3.6485 

12.04788 


KV 

9.4398 

72.9700 

240.9576 


In  a similar  manner,  the  steering  gain,  Ksteer»  is  decreased  as 
a function  of  thrust  acceleration  during  tailoff  in  order  to  prevent 
the  decreasing  torque  capability  of  the  engine  from  forcing  the  actuator 
into  rate  or  position  limits.  After  the  thrust  acceleration  reaches 
its  peak  value,  Ksteer  is  computed  from 
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steer 


steer 


- K 


change  (At  „ ~ 


init  - - 'max 


Kchange  is  defined  as 


change 


‘8teerinit  - Ksteerfinal 


max 


Nominally,  K is  set  to  0.01.  For  the  different  values  of 

steerfinal 

steering  gain  corresponding  to  the  three  autopilot  configurations, 
Kchange  *8 


V 

FCS  change 

2 0.0006 

5 0.00158 

8 0.00256 

2.7  Terminal  steering 

As  previously  discussed,  the  steering  modes  used  in  this  study  are 
unstable  at  the  end  of  the  burn,  the  perturbation  sensitivities  being 
l/Tgo  for  Vg  steering,  and  4/TgQ  for  fuel  depletion  steering.  To  avoid 
this  problem,  the  steering  signal  must  be  modified  as  the  time-to-go 
approaches  zero.  Two  approaches  were  considered  for  this  study: 

(1)  Zero  the  rate  command  vector,  *^CMD* 

(2)  Freeze  the  commanded  thrust  vector,  USF. 

The  latter  method  has  the  slight  advantage  of  having  the  final  commanded 
thrust  direction  oriented  in  the  direction  of  the  last  computed  thrust 

vector  command.  Since  steering  results  in  an  almost  constant  attitude 

maneuver  near  the  end  of  the  burn,  the  advantage  of  either  mode  is  prob- 
ably negligible  for  this  mode.  For  fuel  depletion  steering,  zeroing  the 
commanded  rate  may  contribute  to  the  terminal  stability,  but  if  the  fuel 
depletion  maneuver  continues  to  very  near  the  end  of  the  burn,  this  in- 
crease in  stability  may  be  at  the  expense  of  slight  increases  in  the  Vg 
residuals.  The  guidance  commands  are  not  used  in  the  last  second  of  the 
burn,  but  the  AV  vector  and  the  rate  command  vector  are  continually  com- 
puted, and  this  results  in  a slightly  more  accurate  steering  signal  than 
simply  commanding  a zero  rate.  Therefore,  the  approach  we  will  use  in 
the  study  will  be  to  freeze  the  commanded  thrust  vector,  USF. 
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GUIDANCE  METHODS 


3.1  Position  Offset  Guidance 

As  mentioned  in  Section  1.2,  position  offset  guidance  can  be  used 
during  the  powered  flight  portion  of  the  trajectory  to  correct  for 
inaccuracies  in  the  calculation  of  velocity-to-be-gained  caused  by 
assuming  that  the  velocity  change  is  made  impulsively.  The  development 
of  the  position  offset  guidance  method  in  this  section  is  taken  from 
"A  New  Approach  to  Lambert  Guidance,"  by  Timothy  J.  Brand  of  CSDL. 
(Reference  1) . 

The  traditional  guidance  method  involves  the  pre-maneuver  calcula- 
tion of  a target  offset.  This  calculation  requires  accurate  prediction 
of  the  initial  state  and  the  path  to  be  taken  by  the  vehicle  to  the  tar- 
get. In  the  traditional  method,  an  impulsive  velocity  change  is  used  to 
approximate  the  thrust  phase,  then  the  velocity  required  to  coast  from 
the  initial  position  to  the  terminal  position  in  a specified  time  is 
determined  using  Lamberts'  routine.  Direct  numerical  integration  of 
the  equations  of  motion,  accounting  for  gravity  perturbations,  is  then 
used  to  extrapolate  this  required  velocity  and  initial  position  along 
the  trajectory.  The  first  estimate  of  the  target  offset  is  taken  as 
the  negative  of  the  resultant  miss.  Further  iterations  using  the  off- 
set target  derived  from  the  previous  iteration  as  the  desired  final 
radius  vector  for  the  Lambert  routine  are  generally  required  to  accu- 
rately determine  the  offset  target. 

Unless  the  thrust  maneuver  is  very  short,  the  offset  target 
determined  by  assuming  an  impulsive  velocity  change  can  cause  signif- 
icant terminal  errors  when  used  with  the  Lambert  routine.  This  is 
due  to  the  fact  that  the  compensation  for  gravity  perturbations  is 
based  on  prediction  of  those  effects  over  both  the  thrust  and  the 
subsequent  coasting  periods.  Due  to  the  nonzero  length  of  the  thrust 
phase,  the  actual  trajectory  does  not  follow  the  path  predicted  by 
the  impulsive  velocity  change  assumption,  but  rather  a neighboring  path. 
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The  difference  in  the  perturbing  gravitational  acceleration  between  the 
two  paths  accumulates  over  the  entire  trajectory,  causing  a miss 
at  the  target. 

An  improvement,  then,  to  the  traditional  method  would  be  to  deter- 
mine a coasting  trajectory  that  is  coincident  with  the  actual  trajectory 
at  thrust  cutoff.  In  this  case,  the  difference  in  gravitational  accel- 
erations only  affects  the  trajectory  during  the  thrust  phase,  and  as 
shown  by  Brand  in  Reference  1,  can  usually  be  neglected. 

This  improvement  is  the  position  offset  method,  where  the  initial 
position  is  offset  to  lie  on  a coasting  trajectory  such  that  at  thrust 
cutoff,  the  thrusting  path  and  the  coasting  path  are  coincident  as  shown 
in  Figure  33.  This  pseudo-initial  position  is  then  used  in  the  Lambert 


Figure  33.  Geometrical  relationships  of  perturbation  analysis. 

routine  to  accurately  determine  the  velocity-to-be-gained.  An  important 
advantage  of  this  method  is  that  it  simplifies  the  computation  of  Vg 
between  guidance  cycles.  On  a coasting  trajectory  we  define 


= V ' - V 
req 


Since  V ' = g (r ' ) , we  can  write 


(22) 


where  Ag  is  the  difference  in  the  gravitational  acceleration  between  the 
coasting  trajectories,  powered  trajectories,  and  at  is  the  thrust  accel- 
eration. Brand  gives  numerical  justification  for  neglecting  this  term 
in  Reference  (1),  so  we  have  that 

£ “Aip  (23) 


This  equation  indicates  that  Vg  may  be  accurately  updated  between 
computations  of  Vr  ' by  decrementing  at  the  steering  cycle  rate  by 
the  accelerometer-sensed  AV  over  the  steering  period  Tg.  since  Ag  -*•  0 
as  the  coasting  and  powered  trajectories  become  coincident  at  thrust 
cutoff,  this  computation  becomes  more  accurate  near  the  end  of  the  burn 
when  it  is  important  to  have  an  accurate  Vg  for  thrust  cutoff  calculation. 
The  position  and  velocity  at  thrust  cutoff,  ?co  and  VCQ,  can  be  expressed 
in  terms  of  the  initial  state  rQand  VQ  by  the  following  equations 


V = v + 
CO  o 


/co  _ 

(A.J,  + g(r)  ]dt 
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(24) 


r_„  = r 

co 


fCO  - 

o+J  v 


dt 


(25) 


where  is  the  thrust  acceleration  and  g(r)  is  the  gravitational 
acceleration. 

The  state  at  thrust  termination  can  also  be  expressed  in  terms  of 
the  coasting  trajectory  initial  conditions,  r * and  V0',  as 


co 


J CO 


g(r' )dt 


(26) 


co 


r • 
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V'dt 


(27) 
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where  g ' (r)  is  the  gravitational  acceleration  along  the  coasting  path. 
Then,  the  difference  in  initial  velocity  on  the  coasting  and  thrusting 
trajectories  is 


/CO 

iSj  + 9(r)  - g (r ) ]dt 
o 


The  expression  for  the  thrust  acceleration  of  a vehicle  with 
constant  thrust  T and  mass  flow  rate  in  is 


\ - — ~~  h 

Mq  - mt 


where  i^  is  a unit  vector  in  the  direction  of  the  thrust  vector  and  HQ 
is  the  initial  mass.  From  Equations  24  and  26,  the  velocity  difference 
AV  between  the  velocity  on  the  coasting  trajectory  V*  and  the  powered 
trajectory  ^ can  be  written 


V'  - V 


-f-l± L 

•T  M_  - ml 


where  A7q  is  the  initial  velocity  difference  and  the  gravity  difference 
is  ignored. 


Evaluating  this  integral  gives 


**  = o - Vex  *A  1o* 


1 - t 

T 

. - 


where  t * M0/m  and  the  exhaust  velocity 


Vex  = T/m 


Since  AV  * 0 at  the  cutoff  time 


Avo  = ^ (i  - -v2) 
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We  can  solve  this  equation  for  the  cutoff  time  and  derive 


-AV  /V 

fcco=  - e ° 


Prom  Equations  25  and  27,  the  expression  for  the  initial  position 
offset  is 


rco  _ _ 5co 

ro  = " J (V*  - V)dt  = -J  AV  dt 

O n 


Substituting  for  AV  , we  obtain 


- -rv 


Avo  - vex  H log  (1  - t)  |at 
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Integrating  this  expression 


-ATolco  - T 


vex  H (l  * T2)  **  f1  - -P)  * T2 
- - 


Replacing  Vi.  log  (1  - -SS)  by  AV  , Equation  31  reduces  to 

“X  A TO 


-T  AV  - t Vov  i, 
o co  ex  A 


which  can  be  used  to  determine  the  initial  position  offset  after  the 
cutoff  time,  t , is  computed  from  Equation  30. 

For  fuel  depletion  steering,  the  vehicle  thrust  acceleration  is  not 
pointed  in  the  direction  of  Vg,  but  instead  is  rotated  from  the  Vg 
vector  by  the  amount  of  the  fuel  depletion  angle.  For  this  case 
with  constant  thrust  and  linear  mass  loss,  it  is  shown  in  Reference  3 
that  the  position  offset  equation  is 


Ar  = 


= -K  rv  V + V 2 [cos  0 - v_/v  l B/el 
po|_  cap  g cap  L g'  capj  J 
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where  B is  a unit  vector  normal  to  V and  and  K are  functions  of 

g cap  po 

the  elapsed  time,  t (measured  from  the  time  of  the  initial  state  con- 
ditions) defined  by  the  expressions 


V = I = du 

C3P  U <Mo  - 


(u  - t) 


- p=t 


(Mc  - my) 


Since  the  fuel  depletion  angle  0=0  for  steering,  the  position  offset 
equation  can  be  written  as 


Ar  = -K V V„ 

po  cap  g 


Equations  33  and  35  are  implemented  in  the  guidance  routines  of  the 
simulation  programs  OFFSETT  and  STEER.  For  assumed  constant  thrust,  T, 
and  constant  mass  flow  rate  of  fuel,  m,  the  integrals  in  equation  34  can 
be  solved  to  yield  the  following  solution  for  K 


M - mt 

(t  - t)  log  — + ( 

M - mt 

L o CO  J 

(M  — mt  \ " 
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+ (t  - fcco> 


T o 

where,  as  previously  defined,  V = — and  x = — 

ex  m m 


For  in-flight  mechanization  (as  represented  by  the  STEER  program) , 
the  guidance  computer  solves  for  the  required  velocity  and  the  position 
offset  using  these  equations  at  each  guidance  cycle  time  (every  1.35 
seconds) , and  the  velocity-to-be-gained  is  updated  every  steering  cycle 
(0.45  second)  by  subtracting  from  the  previous  Vg  the  change  in  velocity 
over  the  steering  period.  The  velocity  is  updated  every  autopilot  cycle 
by  accumulating  the  accelerometer  sensed  velocity  increments. 


• * * * 
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3.2  Reentry  Angle  (y)  and  Time-On-Target  (TOT)  Control 

The  fuel  depletion  steering  method  coupled  with  position  offset 

guidance  provides  a convenient  technique  for  simultaneously  controlling 

the  reentry  angle,  y , and  the  Time-On-Target  (TOT) . By  modifying  the 

ratio  V /V  , the  fuel  depletion  angle  may  be  modulated  to  introduce 
g cap 

a pitch  profile  that  modifies  the  trajectory  of  the  missile  such  that 
Y and  TOT  are  controlled  simultaneously. 

Time-On-Target  is  readily  controlled  since  time  of  flight  is  a 
specification  to  the  Lambert  guidance  routine,  i.e.,  the  routine  com- 
putes a trajectory  that  satisfies  the  time  of  flight  and  hence  the 
Time-On-Target  constraint. 

For  reentry  angle  control,  the  guidance  system  first  computes  the 
reentry  angle  as  the  flight  path  angle  at  an  altitude  of  300,000  feet  on 
a trajectory  that  satisfies  the  Time-On-Target  constraint.  If  this  angle 
differs  from  the  desired  reentry  angle,  the  guidance  program  incrementally 
changes  Vcap  to  modify  the  ratio  Vg/Vcap  an(*  thence  the  fuel  depletion 
angle.  This  resulting  change  in  the  pitch  profile  modifies  the  orbit.  A 
new  position  offset  is  then  computed  since  the  offset  is  a function  of 
both  V and  the  depletion  angle,  and  the  Lambert  routine  computes  a 

CAp 

trajectory  that  satisfies  the  time  of  flight  constraint.  A new  reentry 
angle  is  then  computed.  This  process  continues  iteratively  until  the 
reentry  angle  error  is  within  a desired  tolerance. 

3.3  Perturbation  Sensitivity  of  Fuel  Depletion  Guidance 

We  are  interested  in  determing  the  sensitivity  of  the  fuel  depletion 
angle  to  small  deviations  in  the  velocity-to-be-gained  caused  by  guidance 
inaccuracies  and  autopilot/steering  errors.  Obviously,  if  the  computation 
of  Vg  is  exact  and  the  vehicle  is  steered  precisely  along  the  fuel  deple- 
tion arc,  there  will  be  no  stability  problems.  However,  since  there  will 
be  steering  errors,  we  need  to  determine  how  this  sensitivity  varies  as  a 
function  of  the  maneuver  parameters.  We  have  seen  that  in  velocity-to-be- 
gained  steering,  the  sensitivity  of  the  steering  loop  to  steering  errors 

increases  as  1/T  _ and  the  commanded  thrust  vector  USF  must  be  "frozen" 
go 

near  the  end  of  the  burn  to  prevent  instability.  We  expect  to  find  a 
similar  problem  in  fuel  depletion  steering;  as  we  shall  demonstrate,  the 
sensitivity  of  the  steering  loop  in  this  mode  is  actually  4/TgQ 

This  indicates  that  steering  instability  occurs  in  fuel  depletion 
steering  at  a value  of  Tgo  four  times  as  great  as  in  velocity-to-be-gained 
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steering.  This  variable  sensitivity  is  represented  in  Figure  5 in 
Section  2.3  as  the  y/Tgo  term  in  the  steering  loop  model. 

To  analytically  determine  this  sensitivity,  we  introduce  a small 
perturbation,  Verr/  in  the  velocity-to-be-gained  in  a direction  normal 
to  the  fuel  depletion  arc  and  the  unperturbed  thrust  vector  command 
USF^.  Referring  to  Figure  34,  we  call  the  perturbed  vector  V^, 

a0t  = 02-01  -0fl 


'U?F  2 . uTF, 


Figure  34.  Geometrical  relationships  cf  perturbation  analysis. 

and  the  vector  sum  of  Verr  and  Vgl  is  designated  Vg2-  We  assume  that 
this  velocity  error  is  the  result  of  accumulated  steering  errors  such 
that  its  magnitude,  verr»  is  given  by  the  integral 


Verr  -/ 


ATAdt 


where  A is  the  agnle  between  AV  and  the  commanded  thrust  vector  as 
defined  in  Appendix  A3.  If  we  assume  that  the  thrust  acceleration,  A^,, 


is  constant,  we  can  write 


T/M  / Adt 


We  can  write  the  remaining  impulse,  at  this  point  again  assuming 

cop 

constant  acceleration  as 

a A_  T - T/M  T 

cap  t go  ' go 


So  the  ratio  of  V„_/V  is  now 
err  cap 

V /V 
err'  cap 
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From  the  geometry  of  Figure  34,  we  see  that 


V - Verr2  + Vgl2  * 2Vgl  Verr  cos  (90  + V 


writing  Vg22  = (Vgl  + 6Vg) 2 and  cos  (90  + e1)  = - sin  6^  we  now  have 


+ 2v  . 6v  + sy?  = v 2 + 
gl  g s g err 
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+ 2Vgl  Verr  sin  01 


or 


1 Verr 

6Vg  = i -vHr  + Verr  sin  81 
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Also  from  geometry,  we  have  that 


(38) 
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8^  cos  61  - sin  6^ 
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6Vg  Vcap  (61  cos  61  ' Sln  61) 


Expanding  cos  6 and  sin  6 in  power  series  around  6=0  and  subtracting 
gives 

.2 
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and 
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3 6Vq 
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(40) 


Equation  4 gives  the  change  in  6 due  to  a small  variation  in  the  mag- 
nitude of  the  velocity- to-be-gained.  Substituting  for  6Vq  gives 
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Examining  the  two  terms  inside  the  bracket  discloses  that,  while 
both  terms  have  singularities  at  8^  = 0,  the  first  term  is  much  smaller 
than  the  second  term  during  a nominal  burn.  From  a typical  fuel  deple- 
tion simulation  run  (OFFSET),  we  have  that  at  T^q  = 1 second,  V ^ = 

250  ft/sec  and  0^  = 0.113  radian.  For  these  values,  we  see  that  for 

Verr  = 5 ft/sec 


err 
"gl  01 


0.265 
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2.994 


Therefore,  we  neglect  the  contribution  of  the  first  term  and  we  write 
the  expression  for  the  change  in  6 due  to  a variation  in  the  magnitude 
of  the  velocity-to-be  gained  as 
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To  determine  the  total  change  in  the  fuel  depletion  angle  relative 
to  the  unperturbed  vector,  we  must  subtract  the  angle  of  rotation 
0^  since  02  is  computed  relative  to  Vg2  rather  than  Vgl.  This  gives 


total 


V /sin  6,  \ V __ 
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Near  the  end  of  the  burn,  we  can  determine  an  approximate  expression  by 
assuming  that 


sin  6. 


1 ; cos  8^  = 1 


V _ = V 
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We  can  then  write  Equation  43  as 
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We  have  already  shown  that 


Vcap  Tgo  J 


So  now  we  have 
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We  now  wish  to  compare  the  sensitivity  determined  from  the 
approximate  analytical  expression  of  Equation  43  with  the  actual  sen- 
sitivity. To  do  this,  we  introduce  a perturbation  subroutine,  PERT, 
in  the  guidance  program,  OFFSETT.  After  the  nominal  parameters  (6, 
Vcap,  USF^)  are  computed  at  a particular  time,  we  perturb  the 
vector  by  adding  to  v ^ a small  perturbation  vector,  Verr,  and  use 
this  new  Vg2  to  recompute  a new  fuel  depletion  angle  @2  and  a new 
commanded  thrust  vector,  USF2  as  shown  in  Figure  34.  Since  we  are 
computing  these  perturbation  parameters  at  the  same  instant  of  time 
as  the  unperturbed  parameters,  Vcap  remains  constant.  Note  that  the 
magnitude  of  Vg2  is  equal  to  Vgl  plus  a small  increment,  6Vg  and  Vg2 
is  rotated  from  V ^ by  a small  angle  0^.  The  actual  parameters  are 
computed  by  the  PERT  routine  as  follows 


A0  total  = arcsin  jUSFj^  * USF2| 


0 Total 
V 


To  compare  the  actual  sensitivity  to  the  analytical  sensitivity, 
we  first  compare  the  actual  rotation  of  V^2  with  respect  to  Vgi  to  the 
approximation  given  by  Equation  39.  The  actual  rotation  is  computed  by 


arcsin  |Unit(Vgl)  x Unit(Vg2) 


and  the  approximate  rotation  is  again 


(39) 

Repeated 
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The  difference  (normalized  to  Verr)  of  these  two  terms  we  will  call 
ERRl  and 


ERRl  = 


The  assumptions  made  in  this  approximation  are  that  for  8 <<  1, 

sin  8_  = 8_  and  that  cos  8,  = 1. 
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We  now  look  at  Equation  42.  We  have  assumed  that  sin  8^/8^  = 1 

and  that  3/2  <verr/vgi0>  can  be  neglected.  Therefore,  the  approximate 

sensitivity  contribution  of  this  term  is  3 V /v  . The  actual 

err  cep 

contribution  of  this  term  to  the  sensitivity  is  computed  in  PERT  as 
the  difference  between  the  fuel  depletion  angles  computed  from  the 


perturbed  and  unperturbed  Vg  vectors  and  the  constant  vcap. 
then  defined  as 


ERR3  is 


ERR  3 = 


( 6 2 - 0!> 


The  total  sensitivity  of  the  fuel  depletion  angle  to  Vg  perturbations 
is  the  sum  of  ERRl  and  ERR3.  We  therefore  write 


ERR4  * - y— — — (82  - e - 8 ) 

cap  err  y 


We  can  compute  the  sensitivity  partial  by  normalizing  the  total 
Change  in  0 to  Verr  so  we  have 


DTH/DV  = (02  - 81  - 0 )/V( 


g"  err 


The  values  of  these  four  parameters  are  shown  with  the  unperturbed  values 
of  Vg,  Vcap,  and  6 for  a nominal  fuel  depletion  burn  in  Table  13  and 
for  a fuel  depletion  burn  with  reentry  control  in  Table  14.  Note  that 
the  sensitivities  computed  with  and  without  reentry  control  are  es- 
sentially the  same. 
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Table  13.  Comparison  of  analytical  and  computed  sensitivities  of  fuel  depletion 
angle  to  variations  in  the  velocity-to-be-gained  (no  reentry  control) 


SECTION  4 


ATTITUDE  MEASUREMENT 


4.1  Attitude  Measurement  System 

As  previously  discussed,  the  vehicle  uses  the  Advanced  Inertial 
Reference  System  (AIRS)  to  measure  the  attitude  of  the  vehicle  with 
respect  to  an  inertial  reference.  The  AIRS  platform,  shown  in 
Figure  35,  consists  of  a hydraulically  floated  inner  sphere  inside 
an  outer  spherical  case.  The  inner  sphere  contains  the  gyroscopes 
that  are  used  to  maintain  the  orientation  of  the  inner  ball  relative 
to  an  inertial  reference  and  the  accelerometers  which  measure  vehicle 
inertial  accelerations  and  compute  the  velocity  increments.  The  inner 
sphere  is  maintained  at  an  initial  alignment  relative  to  inertial  space 
by  a system  of  hydraulic  jets  that  rotate  the  inner  ball  in  response 
to  signals  from  the  gyros.  The  inertial  ball  has  three  printed  circuit 
resolver  "driver"  bands  which  are  mounted  as  three  orthogonal  great 
circles  on  the  outer  surface  of  the  ball.  The  platform  attitude  is 
measured  in  terms  of  the  intersection  pointB  of  the  driver  bands  with 
a "receiver"  band,  also  mounted  as  a great  circle  on  the  inner  surface 
of  the  case.  The  positions  of  these  intersections  are  defined  in  terms 
of  the  angles  x^>  X2'  X3  measured  along  the  receiver  band  and  in  terms 
of  the  angles  <t>1#  $2,  <t> 3 measured  along  the  driver  bands. 

As  illustrated  in  Figure  36,  the  three  receiver  band  angles  are  all 
measured  from  the  intersection  with  the  receiver  band  of  the  vector  Q^, 
in  the  plane  of  the  receiver  band,  in  a direction  that  would  cause  a 
right-handed  screw  to  advance  along  the  vector  g3  (which  is  perpendi- 
cular to  the  plane  of  the  receiver  band) . The  terminal  points  of  the 
receiver  band  angles  are  defined  as  follows  in  terms  of  the  receiver- 
driver  band  intersections  and  the  triad  of  vectors  P2,  £3,  P^r 
perpendicular  to  driver  bands  1,  2,  3,  respectively:  the  terminal 
point  of  Xi  (where  i - 1,  2,  or  3)  occurs  where  the  receiver  band 
crosses  driver  band  i going  from  the  positive  to  the  negative  side  of 
band  i,  as  defined  by  its  vector 
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The  three  driver  band  angles  4’j_*  ^ » are  measured  from  the 

intersections  with  the  ball  bands  of  the  vectors  £2,  £3,  P^ , re- 
spectively, in  directions  that  would  advance  right-handed  screws  along 
El'  ?-2’  £3'  respectively.  The  terminal  points  of  the  driver  band 
angles  are  where  the  driver  bands  cross  the  receiver  band  from  the 
positive  to  the  negative  side  as  defined  by  its  vector  Q3> 

It  is  important  to  note  that  the  point  of  intersection  of  any 
driver  band  with  the  receiver  band  is  accurately  measurable  only  if  the 
angle  of  intersection  between  these  bands  is  greater  than  45  degrees. 
Fortunately,  at  least  two  of  the  driver  bands  fulfill  this  intersection 
requirement  for  any  ball-case  orientation,  and  only  the  intersections  of 
these  two  driver  bands  with  the  receiver  band  need  be  measured  to  determine 
the  orientation.  The  AIRS  programs  used  in  this  study  are  capable  of 
selecting  those  driver  bands  that  have  the  required  ^45°  intersection 
with  the  receiver  band  and  are  capable  of  determining  the  ball-case 
orientation  from  any  selected  pair  of  driver  band  angles  (e.g.,  4^,  4^) 
in  combination  with  the  corresponding  receiver  band  angles  (e.g.,  Xj)* 

4.2  Attitude  Data  Processing 

The  incremental  updating  of  attitude  errors  following  their 
initialization  at  the  beginning  of  the  Stage  3 burn  is  simple  and  con- 
venient for  an  AIRS  based  autopilot.  This  incremental  approach,  illus- 
trated in  Figure  37,  requires  the  computation  of  an  incremental  trans- 
formation matrix  B which  transforms  a set  of  two  driver  band  angle 
increments  { A<4> ^ , A4>j)  and  one  receiver  band  angle  increment  (Ax^)  into 
the  corresponding  increments  in  body  angles  (AB^,  A&2,  ^e3). 


PITCH  AND  YAW 
f ATTITUDE  ERRORS 

Ezn  ) TO  PITCH  AND 
c ( YAW  AUTOPILOTS 


JDAP  ^BANDAN 

\DAP  A*  VECTOR  I 

EVERY  T, 

lAXiJ  1 


ANGLE  INCREMENT 
)R  DETERMINED 
y Tqap  wcondt 


INCREMENTAL  TRANSFORMATION  MATRIX 


■ AUTOPILOT  SAMPLING  PERIOD 


Figure  37.  Incremental  updating  of  attitude  errors. 
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w = 1/sin  (Xi  - Xj) 


Ideally,  the  B-matrix  should  be  updated  every  autopilot  sampling 
period  prior  to  the  computation  of  each  new  set  of  body-angle  incre- 
ments. However,  the  computation  of  the  trigonometric  functions  of  the 
band  angles  which  make  up  the  elements  of  B is  such  a time-consuming 
process  that  it  is  desirable  to  consider  ways  of  reducing  the  com- 
putation time  by  approximating  the  elements  of  B and/or  updating  these 
elements  less  often  than  every  autopilot  cycle. 
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AIRS BAB  Method 


Various  schemes  for  approximating  the  B-matrix  and  updating  it 
less  frequently  were  studied  at  CSDL  in  1973  and  1974  as  part  of  the 
development  of  the  .single-rotation-axis  (SRA)  autopilot.  This  auto- 

i 

pilot,  which  was  developed  for  large  attitude  maneuvers  of  a space 
vehicle  under  the  control  of  reaction  jets,  takes  advantage  of  a 
i theorem  of  Euler,  according  to  which  the  attitude  of  a body  may  be 

changed  from  one  orientation  to  any  other  orientation  by  rotating 
i the  body  about  an  axis  which  is  fixed  to  the  vehicle  and  stationary 

1 in  inertial  space.  The  SRA  autopilot  study  showed  that  if  the  values 

of  the  band  angles  used  in  the  B-matrix  updates  are  extrapolated  every 
0.45  seconds  into  the  middle  of  the  next  0.45-second  time  interval,  the 
B-matrix  based  on  these  extrapolated  angles  can  be  employed  over  the 
entire  0.45-second  interval,  with  negligible  degradation  in  SRA  auto- 
pilot performance  compared  to  updating  B precisely  every  0.03  seconds. 
This  method  of  updating  the  B-matrix  has  been  termed  "AIRSBAE"  (“AIRS 
| band  angle  extrapolation”). 

The  AIRSBAE  method  is  illustrated  in  Figure  38,  where  the  matrix 
! updating  interval  is  designated  as  Ts>  Here,  it  should  be  noted  that 


Extrapolated  value  of  band  angle 


(n  • DT,  (n-l)T,  nT$  (n  + l)Tt  <"  + 1>Ts 

2 * . 2 * 

(now) 


Figure  X.  Pictorial  representation  of  AIRSBAE. 
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the  same  symbol  (Ts)  is  used  for  both  the  steering-loop  sampling  period 
and  the  B-matrix  updating  interval  because  the  studies  carried  out  in 
this  thesis  and  in  the  SRA  autopilot  investigations  have  shown  that  the 
choice  of  Tg  • 0.45  seconds  for  both  the  sampling  period  and  updating 
interval  is  reasonable  for  the  dynamic  requirements  considered.  This 
period  amounts  to  fifteen  of  the  30  ms  autopilot  sampling  periods  which 
are  used  in  this  thesis  (as  well  as  in  the  earlier  SRA  autopilot  studies) . 
However,  the  investigations  of  B-matrix  updating  techniques  described 
below  consider  also  other  values  of  the  B-matrix  updating  interval  of 
2,  30,  and  45  autopilot  cycles,  which  may  not  be  the  same  as  the  optimum 
sampling  period  for  steering- loop  computations . Thus,  the  symbol  T^ 
can  represent  a different  time  interval  in  the  case  of  the  B-matrix 
updating  than  in  the  case  of  steering. 

Referring  to  Figure  38,  it  is  seen  that  the  AIRSBAE  method  pre- 
dicts the  band  angles  halfway  through  the  next  Tg  interval  by  taking 
half  the  difference  between  the  present  value  of  each  band  angle  and 
the  value  measured  Tg  seconds  previously,  and  adding  that  difference 
to  the  present  value  of  the  band  angle. 

The  dead  zones  and  allowed  terminal  errors  in  both  angular  velocity 
and  attitude  in  the  SRA  autopilot  provided  some  tolerance  for  the  errors 
resulting  from  the  infrequent  updating  of  the  B-matrix  in  the  AIRSBAE 
approach.  However,  the  dynamic  requirements  of  other  autopilots, 
such  as  the  boost  autopilot  being  considered  in  this  thesis,  may  not 
permit  the  errors  in  attitude  and  estimated  angular  rate  which  result 
from  the  AIRSBAE  approach. 

One  aspect  of  particular  concern  in  applying  the  AIRSBAE  method 
to  the  third  stage  boost  autopilot  with  fuel  depletion  and  gamma -time-on- 
target  guidance  is  the  relatively  large  step  change  in  the  B-matrix 
occurring  every  0.45  seconds  as  the  result  of  the  large  angular  rates 
associated  with  these  guidance  methods.  This  step  change  in  B can  pro- 
duce autopilot  transients  which  might  not  be  considered  acceptable  in 
these  boost  applications. 

AIRSME4  Method 

The  AIRSBAE  method  was  considered  along  with  other  B-matrix  up- 
dating approaches  in  terms  of  the  Stage  3 boost  autopilot  and  the 
dynamic  requirements  imposed  by  fuel  depletion  guidance  and  gamma-time- 
on-target  guidance.  The  most  promising  alternative  to  the  AIRSBAE 
technique  has  been  termed  "AIRSME4"  ("AIRS  matrix  extraploation  4"). 
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This  alternative,  which  is  illustrated  in  Figure  39,  extrapolates  the 
band  angles  to  the  end  of  the  next  T&  interval  (instead  of  to  the  middle. 


actual  band  angle 


extrapolated 
band  angle 


B-Matrlx 


BAND 

ANGLE 


incramentally 
••ch  DAP 
cyde 


N(nTt) 
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N (n  + 1)T, 


(n  - 2)T$ 


<n-1>Tt 


(n  ♦ 1)Tg 


at(n+1)Tt."o!d"N(n+1)T| 
ax  trap 


N (nTf) 
extrap 


Figure  39.  Pictorial  representation  of  A IRS HE 4 . 

aa  in  AIRSBAE) , and  computes  a predicted  B-matrix  from  these  extrapolated 
angles.  The  predicted  B-matrix  from  the  previous  T^-interval  is  then 
subtracted  from  the  new  predicted  matrix  and  the  resulting  matrix  dif- 
ference divided  by  the  number  of  autopilot  cycles  in  T(  to  obtain  an 
incremental  B-matrix,  AB.  The  B-matrix  is  then  updated  every  autopilot 
cycle  within  T(  by  adding  AB  to  the  B-matrix  determined  in  the  pre- 
vious cycle,  starting  with  the  predicted  B-matrix  from  the  previous  T# 
interval. 

The  band  angle  extrapolation  approach  used  in  AIRSMF4  is  based  on 
differencing  the  band  angles  over  T^,  as  is  done  in  AIRSBAE,  but  in 
AIRSME4  the  entire  difference  is  added  (rather  than  half  the  difference) 
to  perform  the  extrapolation. 

The  incremental  updating  of  the  B-matrix  every  autopilot  cycle  by 
the  AIRSME4  approach  results  in  a smoother  operation  of  the  autopilot 


tfWSfc-  A.*- 


K * • 4 


than  with  AIRSBAE,  and  accomplishes  this  improvement  with  only  a small 
increase  in  the  computational  burden. 

Band  Angle  Switchover 

The  switching  of  band  angles,  which  can  occur  of  any  autopilot 
cycle,  is  treated  in  the  following  manner  as  a special  case  in  both 
AIRSBAE  and  AIRSME4 : 

(1)  The  B-matrix  is  recomputed  in  terms  of  the  new  set  of  band 
angles  as  soon  as  the  switchover  occurs. 

(2)  The  switchover  value  of  the  B-matrix  is  retained  for  a 
pre-specif ied  number  of  autopilot  cycles,  uSWitcj,over ' 
which  is  equal  to  5 in  all  the  simulation  results  pre- 
sented in  this  thesis. 

(3)  When  uswitch0ver  autopilot  cycles,  have  elapsed,  the 
band  angles  are  again  extrapolated  (to  the  middle  of 
Tg  in  AIRSBAE  and  to  the  end  of  T^  in  AIRSME4)  and 

the  standard  B-matrix  updating  approach  is  resumed.  In 
the  case  of  AIRSME4  the  value  of  B needed  at  the  be- 
ginning of  the  new  T interval  is  computed  from  current 

8 

band  angles. 


Computation  Requirements 

A comparison  of  the  computation  times  required  by  AIRSBAE  and 
AIRSME4  for  an  updating  time  interval  of  Tg  » 0.45  sec  is  presented  in 
Table  15.  Here,  it  is  seen  that  AIRSBAE  requires  4.65  ms  per  0.45  sec, 
while  AIRSME4  requires  5.38  ms  per  0.45  sec.  Thus,  the  difference  in 
computation  times  of  the  two  updating  methods  is  an  almost  negligible 
value  of  0.73  ms  per  0.45  sec  (which  is  less  than  0.05  ms  per  30  ms 
autopilot  cycle).  However,  it  should  be  pointed  out  that  both  B-matrix 
techniques  are  very  substantial  improvements  over  updating  the  B-matrix 
from  its  trigonometric  functions  every  30  ms  autopilot  cycle,  which 
would  require  about  36  ms  per  0.45  sec.  The  reduction  from  36  ms  per 
0.45  sec  to  an  approximate  time  of  5 ms  for  AIRSBAE  and  AIRSME4  amounts 
to  a reduction  from  an  8%  computation  load  to  a 1.1%  computation  load 
for  attitude  data  processing. 


4.3  Idealized  AIRS  Error  Study 

The  first  phase  of  the  attitude  measurement  study  evaluated  the 
purely  geometric  errors  resulting  from  AIRSBAE  and  AIRSME4  when  the 
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Table  15.  Comparison  of  computational  requirements  of  AIRSBAE  and 
AIRSME4*  based  on  the  Honeywell  701P  computer. 


r 


Li 


Operations  Common  to  Both  AIRSBAE  and  AIRSME4 


Computation  Time 
s 


1. 

Update  Incremental  Transformation  Matrix 
at  Beginning  of  Tfi 

Units 

(m  sec)  per 

4 sin/cos  § 125  units** 

m 

500 

2.23 

14  Mult./Div  @ 3 units 

a 

42 

0.19 

1 Subtr.  9 1 unit 

a 

1 

0.0045 

2. 

Matrix/Vector  Multiplication  at  Each 
Autopilot  Cycle  (^DAp) 

9 Mult.  # 3 units  x 15 

as 

405 

1.80 

6 Add.  8 1 unit  x is 

a 

90 

0.40 

Totals 

1038 

4.62 

AIRSBAE 

1. 

Common  Operations 

as 

1038 

4.62 

2. 

Extrapolate  Band  Angles  Over  Tg 

6 Add.  9 1 unit 

m 

6 

0.03 

(Division  by  2 is  handled  by  a SHIFT 
operation  which  takes  neglible  time) 

Totals 

1044 

4.65 

AIRSNE4 

1. 

Coeuon  Operations 

m 

1038 

4.62 

2. 

Difference  Incremental  Matrices  at  T( 

9 Subtr.  9 1 unit 

m 

9 

0.04 

3. 

Division  of  Tg  Increments  by  TDAp 

9 Div.  8 3 units 

m 

27 

0.12 

4. 

Addition  of  TDAp  Increments  to  Matrix 
at  Each  Autopilot  Cycle 

9 Add  9 1 unit  x is 

- 

135 

0.60 

Totals 

1209 

5.38 

•Neglecting  Computer  Bookkeeping  Operations,  assuming  T ■ 15  autopilot 

cycles 


•*1  Unit  “ 4.45  usee. 
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» 0.45  sec 


boost  veohicle  is  rotating  at  constant  and  equal  rates  about  the  pitch 
and  yaw  axes,  with  no  motion  about  the  roll  axis.  Two  programs  were 
developed  for  AIRSBAE  and  AIRSME4  which  determine  the  errors  in  com- 
puted body-angle  increments  per  degree  of  separate  pitch,  yaw  and  roll 
rotations  which  occur  when  the  vehicle  rotates  at  specified  rates 
at  various  orientations.  Vehicle  orientations  considered  were  defined 
in  terms  of  roll  Euler  angles  of  0°,  45®,  90®,  135®,  and  180®  pitch 
Euler  angles  which  are  incremented  by  30®  from  -180®  to  -1-180®  and 
yaw  Euler  angles  which  are  incremented  by  30®  from  -90®  to  +90®.  These 
programs  also  search  for  the  worst-case  errors  in  body-angle  increments 
per  degree  rotation  for  each  pair  of  equal  pitch  and  yaw  angular  rates 
assumed.  Rates  were  considered  which  yield  pitch  and  yaw  rotations 
over  Tfl  seconds  of  2.5®,  5®,  10®,  and  20®.  The  programs  compute  the 
body-angle  errors  at  only  those  points  where  there  is  not  a switchover 
of  the  attitude  bands  over  the  Tg  intervals  used  for  computing  these 
errors.  A comparison  of  the  worst-case  body-angle  errors  per  degree 
rotation  about  any  body  axis  is  given  in  Table  16.  This  comparison, 
indicates  consistently  lower  errors  for  AIRSME4. 

4.4  Autopilot  Simulation  Study  Based 
on  the  STAR  Program 

The  next  step  in  the  evaluation  of  the  attitude  measurement  study 
was  to  use  AIRSBAE  (renamed  AIRS)  and  AIRSME4  (renamed  AIRSX)  as  sub- 
routines in  the  STAR  simulation  program.  The  STAR  program  is  basically 
a simulation  of  the  Stage  III  vehicle  operating  in  inertial  space  with 
no  guidance  system.  The  vehicle  is  steered  by  commanding  a pitch  rate 
profile  that  approximates  the  profile  seen  in  fuel  depletion  guidance. 
This  profile  consists  of  a parabolically-increasing  pitch  rate  command 
until  the  pitch  rate  reaches  31.25  deg/sec,  then  a constant  31.25 
deg/sec  rate  is  commanded  for  two  seconds.  At  that  time  a linearly- 
decreasing  rate  signal  is  commanded  for  two  seconds,  followed  by  a 
steady  state  10  deg/sec  rate  signal.  This  commanded  pitch  profile  is 
shown  on  channel  two  of  the  response  curves.  The  displayed  variables 
are  given  below. 

PITCH  RATE  - The  vehicle  pitch  rate  in  deg/sec. 

PITCH  RATE  COMMAND  - The  commanded  pitch  profile 

INDICATED  RATE  - The  vehicle  pitch  rate  as  indicated  by  the 
AIRS  (AIRSBAE)  or  AIRSX  (AIRSME4)  platform. 
In  a NOAIRS  run,  this  measured  rate  is  not 
used  as  the  feedback  signal. 
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ENGINE  DEFLECTION  - The  engine  actuator  pitch  deflection  in 

degrees. 

PITCH  ATTITUDE  - The  vehicle  pitch  attitude  in  degrees. 

REJECTED  BAND  - The  driver  band  that  has  an  intersection  with 
a receiver  band  of  less  than  45  degrees. 

Since  only  two  band  intersections  are  required 
to  completely  determine  the  orientation  of 
the  AIRS  inner  sphere  relative  to  the  case, 
this  band  is  not  needed  and  is  rejected. 

Plots  of  the  above  variables  versus  time  are  presented  in  Figures 
40  through  55  for  various  combinations  of  the  following  conditions: 

(1)  Type  of  feedback: 

(a)  Ideal  feedback  (denoted  as  "NOAIRS")  in  which  the  actual 
body-angle  increments  are  employed  for  feedback. 

(b)  AIRS-based  feedback  in  which  the  feedback  body-angle 
increments  are  computed  either  by  an  "AIRS"  subroutine 
(based  on  the  AIRSBAE  approach)  or  by  an  "AIRSX"  sub- 
routine (based  on  the  AIRSME4  method) . 

(2)  AIRS  errors: 

(a)  If  the  types  of  errors  in  the  AIRS  band-angle  measure- 
ments are  not  specified  in  the  title  of  the  simulation 
run,  these  errors  have  been  assumed  zero. 

(b)  The  measured  angle  of  any  driver  or  receiver  band  con- 
tains a deterministic  error  term  which  is  a function 
of  the  angles  of  both  the  driver  and  receiver  band 
angles  which  define  any  particular  band  intersection. 

The  maximum  magnitude  of  this  error  term  is  less  than  ten 
arc  minutes  although  its  rate  of  change  can  be  as 

large  as  3 degrees  per  degree  of  band  angle  change.  A 
table  lookup  representation  of  the  deterministic  errors 
based  on  laboratory  measurements  is  employed  when  these 
errors  are  simulated. 

(c)  A second  source  of  band  angle  measurement  errors  is  the 
random  electrical  noise.  This  noise  is  represented  in 
simulation  runs  by  a random  number  generator  whose 
standard  deviations  are  pessimistically  assumed  to  be 
180  arcseconds  for  the  driver  bands  and  20  arc  seconds 
for  the  receiver  bands. 
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(3)  Time  interval  between  updates  of  B based  on  extrapolated 
band  angles: 

The  number  of  30  ms  autopilot  cycles,  UMAX,  between 
updates  of  the  B-matrix  based  on  extrapolated  band  angles 
is  considered  to  have  the  following  values  in  the  simulation 
runs:  2,  15  (nominal  case),  30  and  45. 

The  simulation  runs  in  Figures  40  through  55  are  described  in 
terms  of  the  above  conditions  in  Table  17. 

It  is  seen  in  the  figures  of  simulation  runs  that  increasing  UMAX 
past  15  cycles  gives  a significant  decrease  in  the  fidelity  of  the  feed- 
back signal  for  both  AIRS  and  AIRSX  runs.  For  all  UMAX  values  the 
AIRSX  data  is  of  significantly  higher  quality  than  the  AIRS  data.  This 
is  because  the  incremental  matrix  is  updated  at  each  autopilot  cycle 
in  AIRSX,  rather  than  at  intervals  of  UMAX  autopilot  cycles  as  in  AIRS. 
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Table 

17.  STAR 

simulation 

runs. 

Run  No. 

Mode 

Ideal 

Feedback 

UMAX 

AIRS 

Noise 

Deterministic 

Errors 

4/00 

NOAIRS 

yes 

15 

no 

no 

4/40 

NOAIRS 

yes 

15 

no 

no 

4/21 

AIRS 

no 

2 

no 

no 

4/22 

AIRS 

no 

15 

no 

no 

4/23 

AIRS 

no 

30 

no 

no 

4/24 

AIRS 

no 

45 

no 

no 

4/25 

AIRS 

no 

15 

yes 

no 

4/26 

AIRS 

no 

15 

no 

yes 

4/27 

AIRSX 

no 

2 

no 

no 

4/28 

AIRSX 

no 

15 

no 

no 

4/29 

AIRSX 

no 

30 

no 

no 

4/30 

AIRSX 

no 

45 

no 

no 

4/31 

AIRSX 

no 

15 

yes 

no 

4/32 

AIRSX 

no 

15 

no 

yes 

4/33 

AIRS 

no 

15 

yes 

yes 

4/34 

AIRSX 

no 

15 

yes 

yes 
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SECTION  5 


SIMULATION  RESULTS 


5.1  Overview 

This  section  presents  the  results  of  the  computer  simulations 
using  the  OFF SETT  and  STEER  programs.  As  explained  in  Appendix  C, 

OFF SETT  is  the  simulation  of  the  guidance  system  with  a point  mass  vehicle 
model  with  no  autopilot  or  steering  dynamics.  STEER  is  a complete 
system  simulation  that  incorporates  the  guidance  equations  from  the 
OFFSETT  program,  the  steering,  autopilot,  and  actuator  dynamics  from 
the  STAR  program,  and  the  attitude  measurement  effects  from  AIRS. 

The  simulation  initial  conditions  were  the  same  for  both  the 
OFFSETT  and  STEER  simulation  runs.  These  initial  conditions  are  given 
in  Table  18  along  with  the  results  of  the  solution  of  the  Lambert  rou- 
tine, position  offset  computations,  and  the  fuel  depletion  guidance 
equations  from  the  initial  guidance  cycle. 

Table  18.  OFFSETT  and  STEER  simulation  initial  conditions. 

Altitude  Velocity  Range  Angle  Vcap  Time  of  Flight 
(ft)  (ft/sec)  (deg)  (ft/sec)  (sec) 

400,000  18,000  100  8400  2500 

The  inital  solution  to  the  Lambert  routine  with  position  offset 
guidance  gives  the  following  results  for  the  given  initial  con- 
ditions and  fuel  depletion  steering: 


Required 

Velocity 

(ft/sec) 


Velocity-to-be-Gained 
(ft/sec) 


Reentry 

Angle 

(deg) 


Position 

Offset 

(ft) 


24,070 


6,071 


-27.54 


231,456 


Fuel  Depletion  Angle  (deg) 


5.2  OFFSETT  Simulation  Runs 


The  results  of  the  OFFSETT  simulation  are  given  in  Figures  56  through 
60,  and  the  parameters  of  each  of  the  five  runs  are  given  in  Table  19. 
Figure  56  is  an  example  of  steering,  and  Figure  57  shows  fuel  deple- 
tion steering.  Figures  58-1  and  58-2  show  the  effect  of  simultaneous 
control  of  Time-On-Target  and  reentry  angle.  Figures  59-1  and  59-2 
show  the  effect  of  velocity  perturbations  on  fuel  depletion  steering, 
while  Figures  60-1  and  60-2  show  the  effect  of  velocity  perturbations 
on  fuel  depletion  steering  with  reentry  control. 


Table  19.  OFFSETT  simulation  runs. 


Figure/Run 

Steering 

Mode 

Reentry 

Control 

Perturbation 

Analysis 

56 

2/00 

vg 

No 

NO 

57 

2/01 

Fuel 

Depletion 

No 

NO 

58 

2/03 

Fuel 

Depletion 

Yes 

No 

59 

2/20 

Fuel 

Depletion 

No 

Yes 

60 

2/21 

Fuel 

Depletion 

Yes 

Yes 

The  output  parameters  shown  in  these  figures  are  explained  below: 
DEPL  ANGLE  - The  fuel  depletion  angle,  0 

VCAPT  - The  AV  capability 

VG  MAG  - The  magnitude  of  the  velocity-to-be-gained 

RENTRY  ANGLE  - The  reentry  angle,  y 

RENTRY  ERR  - The  reentry  angle  error,  y - Ydeaired 

ANGL  - The  angle  between  the  perturbed  and  unperturbed 

V_  vectors 

g 

THR  ACCEL  - The  thrust  acceleration 

KPOT  - The  position  offset  constant, 
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VGX 

- The  horizontal  component  of  the  Vg  vector  relative 
to  inertial  axes 

VGZ 

- The  vertical  component  of  the  Vg  vector  relative 
to  inertial  axes 

USFX 

- The  horizontal  component  of  the  thrust  vector  com- 
mand relative  to  inertial  axes 

USFZ 

- The  vertical  component  of  the  thrust  vector  com- 
mand relative  to  inertial  axes 

ERR  1,  3,  4, 

- Differences  between  theoretical  and  computed  per- 
turbation sensitivities 

DTHTOT 

- The  theoretically  determined  partial,  30/3V^ 

THPPRT 

- The  actual  value  of  the  partial,  30/3Vg 

5*3  STEER  Simulation  Runs 

The  results  of  the  STEER  simulation  study  are  given  in  Figures 
61-83.  Examples  of  velocity-to-be-gained  steering  are  given  in 
Figures  61-64,  and  Figures  65-83  illustrate  fuel  depletion  steering. 

The  initial  conditions  for  the  STEER  simulations  are  given  in  Table  18, 
and  the  input  parameters  for  each  of  the  STEER  simulation  runs  are 
given  in  Table  20.  In  Table  20,  the  autopilot  and  steering  con- 
figuration for  each  run  are  listed  in  the  second  column  according 
to  the  following  code.  The  first  number  indicates  the  autopilot  cross- 
over frequency  and  the  number  after  the  slash  indicates  the  steering 
model,  e.g.,  the  5 rad/sec  autopilot  (FCS5)  and  s eering  model  01 
(no  integral  feedback)  are  designated  as  5/01.  nn  asterisk  after  the 
autopilot  frequency  means  that  the  "optimized"  autopilot  and  steering 
combination  is  used  in  that  run.  The  fifth  column  indicates  which  (if  any) 
AIRS  method  is  used.  The  code  0 means  that  ideal  measurements  are  used, 
the  code  1 stands  for  AIRS,  and  code  2 stands  for  AIRSX.  The  tailoff 
Column  indicates  whether  or  not  the  thrust  tailoff  model  is  used.  The 
comments  column  are  self-explanatory. 


117 


Sfeta.-,  Miiw, 


ERr> (DEC  r»TCH  WEID/5ECJ  p R9TE  CHDlD/'SEL)  USFZ  VCZ  IFT/SEC)  VCX  tFf/rjECi 


Figure  63-1.  STEER,  Vg  steering  with  AIRSX  and  thrust  tailoff. 
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PITCH  flfiTE  (D/SEC) 


Sff.tR  ERRIDEC  P'TCH  RRTE  lO/SECJ  P WE  CMD10/5EU  USFZ  VCE  (FT /SEC!  VC/  ifT-'Sf; 


Figure  76-2.  STEER,  fuel  depletion  steering,  8 rad/sec  autopilot,  integral 


5 TEER  ERR  (DEC  PITCH  WE  (D/SEO  p RATE  CMD10/5ECJ  USfZ  VG2  If  T /SEC)  VCX  iFT/SFC) 


STEER  ERR  (OEG  PITCH  RPTE  (0/SEC!  P RRTE  CHD (D/SED  USFZ  VGZ  (FT 'SEP  VCX  iFT/SEC) 
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Table  20.  STEER  simulation  runs. 


Figure/Run 

Autopilot/ 

Steering 

Mode 

T 

FR 

AIRS 

Tailoff 

Comments 

61 

1/01 

5/01 

V 

9 

0 

0 

No 

Normalised  Design 

62 

1/02 

5/01 

V 

9 

0 

2 

No 

Normalized  Design 

63 

1/03 

5/01 

V 

9 

0 

2 

Fes 

Normalized  Design 

64 

1/04 

5/01 

V 

9 

5 

2 

Yes 

Terminal  Steering 

65 

1/05 

5/01 

FD 

0 

0 

No 

Normalized  Design 

66 

1/06 

5/02 

FD 

0 

0 

No 

Integral  Feedback 

67 

1/07 

5/02 

FD 

2 

0 

No 

Terminal  Steering 

68 

1/08 

5/02 

FD 

2 

2 

No 

69 

1/09 

5/01 

FD 

2 

2 

Yes 

70 

1/10 

5/01 

FD 

2 

0 

No 

Simplified  0 Computation 

71 

1/11 

5/01 

FD-Reentry  Ctl. 

2 

0 

No 

Gamma/TOT 

72 

1/12 

5/02 

FD-Reentry  Ctl. 

2 

0 

No 

Gamma/TOT 

73 

1/13 

2/01 

FD 

2 

0 

No 

FCS2 

74 

1/14 

2/02 

FD 

2 

0 

No 

FCS2 

75 

1/15 

8/01 

FD 

2 

0 

NO 

FCS8 

76 

1/16 

8/02 

FD 

2 

0 

No 

FCS8 

77 

1/17 

5*/01 

FD 

2 

0 

NO 

Optimized  Design 

78 

1/18 

5*/02 

FD 

2 

0 

No 

Optimized  Design 

79 

1/19 

5*/02 

FD 

2 

2 

No 

Optimized  Design 

80 

1/20 

5*/02 

FD 

2 

2 

Yes 

Optimized  Design 

81 

1/21 

5*/01 

FD 

2 

2 

Yes 

Optimized  Design 

82 

1/22 

8/01 

FD 

2 

1 

No 

FCS8  w/AIRS 

83 

1/23 

8/01 

FD 

2 

2 

NO 

FCS8  w/AIRSX 

The  output  parameters  shown  in  Figures  61-83  and  not  previously 
explained  (Section  5.1)  are  discussed  below. 


P RATE  COMMAND  - 


The  pitch  rate  command  to  the  autopilot  from 
the  steering  loop 


PITCH  RATE  - The  vehicle  pitch  rate 

STEERERR  - The  angular  difference  between  the  vectors 

USF  and  AV 


PITCH  ATT 

ENG  DEFL-P 
ENG  RATE 


The  vehicle  pitch  attitude  measured  from  the 
local  vertical 

The  pitch  deflection  of  the  engine  actuator 

The  pitch  deflection  rate  of  the  engine 
actuator 


Note  that  in  these  figures,  a negative  pitch  rate  command  or  response 
is  actually  nose-up  due  to  the  sign  convention  established  in  Section  1. 
A negative  pitch  response  (up)  corresponds  to  a negative  rotation  (up) 
of  the  engine  nozzle.  Pitch  attitude,  however,  is  defined  in  the 
opposite  (and  more  logical)  sense.  Negative  pitch  rotation  of  the 
vehicle  will  result  in  a positive  increase  in  the  pitch  attitude,  which 
is  measured  from  the  local  vertical.  Note  also,  that  the  fuel  depletion 
angle  is  measured  from  the  Vg  vector  rather  than  the  inertial  axes. 
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SECTION  6 


CONCLUSIONS 


6 . 1 Overview 

The  purpose  of  this  study  was  to  examine  the  design  of  the  auto- 
pilot and  steering  loops  and  the  attitude  processing  requirements  for 
the  third  stage  of  a ballistic  missile.  In  accomplishing  this,  we  have 
discussed  the  design  of  the  steering  loops  and  examined  the  stability 
of  these  loops  as  a function  of  the  changing  vehicle  parameters.  Then, 
we  developed  a simplified  method  for  determining  the  parameters  of  the 
autopilot  filter  for  any  desired  crossover  frequency  by  normalizing 
the  design  to  a convenient  parameter. 

We  have  examined  the  effect  of  the  steering  loop  on  the  system 
response  and  then  used  the  steering  dynamics  to  "optimize"  the  closed- 
loop  response  in  a simple  manner. 

The  guidance  equations  for  position  offset  guidance  for  velocity- 
to-be-gained  and  fuel  depletion  steering  were  then  discussed,  and  a 
perturbation  study  was  performed  to  determine  analytically  and  em- 
pirically the  sensitivity  of  these  methods  to  steering  errors,  and 
evaluate  the  analytical  expression  so  determined.  Next,  several 
attitude  processing  methods  were  developed  and  compared  as  to  their 
effectiveness  in  reducing  the  computation  load  of  the  attitude  s'  stem. 
Then,  two  of  the  methods  were  evaluated  by  using  them  in  a simulation 
of  the  vehicle  in  maneuvers  involving  high  timing  rates. 

Finally,  the  autopilot  and  steering  loops  and  the  "best"  attitude 
data  processing  method  were  combined  with  the  guidance  methods  and 
evaluated  in  a full  system  simulation  of  a ballistic  missile  flight. 

In  this  simulation,  the  effects  of  data  measurement  noise,  thrust 
tailoff,  and  terminal  steering  were  interjected  to  test  the  system 
response  and  performance  compared  to  an  ideal  system.  The  conclusions 
arrived  at  through  this  process  are  presented  in  this  section. 
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From  the  study  done  to  compare  the  worst  case  errors  for  AIRS 
and  AIRSX,  we  see  that  AIRSX  (AIRSME4)  has  a distinct  advantage  over  the 
other  methods  of  processing  the  attitude  data.  The  worst  case  errors 
for  AIRSX  were  significantly  lower  than  those  for  "standard"  AIRS 
(AIRSBAE) . It  we  also  compare  the  computational  requirements  required 
for  each  method,  we  notice  that,  although,  AIRSX  requires  slightly  more 
computation  time  than  AIRS,  the  difference  in  the  computational  require- 
ments is  negligible  compared  to  the  total  time  required.  This  is  part- 
icularly true  if  the  difference  is  compared  to  the  other  computational 
requirements  of  the  guidance  computer.  When  these  two  methods  are  com- 
pared in  the  STAR  simulation,  we  see  that,  again,  the  AIRSX  data  is 
consistently  more  accurate  than  the  AIRS  data,  particularly  at  the 
higher  rotation  rates. 


6 . 3 STEER  Simulation  Studies 


This  section  presents  the  conclusions  drawn  from  the  OFFSETT 
and  STEER  programs.  We  will  discuss  the  results  of  simulation  runs 
that  were  made  with  the  autopilot  designed  using  the  normalized  design 
process,  and  compare  this  system  response  with  that  of  the  system  using 
the  "optimized"  autopilot  and  steering  combination.  Then  these  results 
will  be  compared  to  the  ideal  case,  OFFSETT. 

Normalized  Design  - The  STEER  simulation  has  verified  that  the 
normalized  design  process  may  be  effectively  used  to  determine  the 
autopilot  parameters  for  any  desired  crossover  frequency.  For  each 
of  the  3 designs  (2,  5,  and  8 rad/sec  autopilot  crossover/frequencies), 
the  frequency  response  plots  verify  that  the  compensated  system  has 
the  desired  open  and  closed  loop  response,  and  the  STEER  simulations 
have  indicated  that  these  systems  behave  as  they  should.  Specific 
conclusions  that  may  be  drawn  from  the  simulation  runs  that  involve 
the  normalized  design  autopilots  are  described  below: 

I.  Steering  Model  - For  fuel  depletion  steering,  the  simplest 
steering  model,  i.e.,  the  proportional  steering  mode,  gave 
the  lowest  residual  velocity-to-be-gained.  Since  this 
model  is  less  accurate  than  it  would  be  if  integral  feedback 
was  added  (as  is  shown  by  the  magnitude  of  the  steering  errors 
for  these  cases) , the  excess  AV  capability  is  wasted  early 
in  the  maneuver,  the  guidance  system  switches  back  to  the  V 


steering  mode  at  an  appreciable  value  of  T 


Hence , more 


time  is  spent  nulling  the  Vg  errors  with  this  steering  model 
than  there  is  with  the  integral  feedback  added. 

II.  Terminal  Steering  - Terminal  steering  is  required  during  the 
last  few  seconds  of  the  burn  when  fuel  depletion  steering  is 
used.  The  instability  is  more  noticeable  when  the  switch 
back  to  Vg  steering  occurs  in  the  last  few  seconds.  Even 
with  terminal  steering  in  the  last  two  seconds,  the  final 
pitch  rates  and  residual  Vg's  are  higher  when  the  fuel 
depletion  maneuver  extends  to  the  end  of  the  burn. 

III.  Simplified  6 Calculation  - The  use  of  the  simplified  expres- 
sion for  computing  6 gave  no  noticeable  degradation  in  re- 
sponse when  used  with  the  simple  steering  mode.  Since  this 
method  always  computes  a smaller  6 than  the  iterative  method, 
particularly  at  the  higher  values  of  6 , when  used  with  a more 
accurate  steering  method  it  could  cause  underwasting  and  the 
fuel  depletion  angle  could  suddenly  increase  at  the  end  of 
the  burn. 

IV.  Attitude  Data  Processing  - The  use  of  AIRSX  or  AIRS  for 
attitude  data  processing  interjects  about  10  deg/sec  of 
actuator  rate  noise  with  5 rad/sec  autopilot  and  about  40 
deg/sec  of  noise  with  the  8 rad/sec  autopilot.  For  the  8 
rad/sec  autopilot,  this  noise  also  causes  significant  resid- 
ual Vg  errors.  Without  AIRS  or  AIRSX,  however,  the  8 rad/sec 
autopilot  gave  the  best  performance  in  fuel  depletion  steering, 
particularly  with  the  simple  steering  model. 

V.  Gamma/TOT  Control  - With  the  5 rad/sec  autopilot  it  was  pos- 
sible to  control  the  reentry  angle  within  the  desired  tol- 
erance (10-6  radians)  from  the  beginning  of  the  burn  until 
t = 49  seconds  for  the  proportional  steering  model,  and  until 
t = 59  seconds  for  the  proportional  plus  integral  model.  For 
the  given  initial  conditions,  however,  an  additional  1600 
fps  of  AV  capability  was  required  to  accomplish  this. 

IV.  Specific  Autopilot  Capabilities  - From  examination  of 

Figures  73  and  74,  it  should  be  clear  that  the  2 rad/sec 
autopilot  is  not  capable  of  accurately  following  the  pitch 
profile  commanded  by  the  guidance  system  during  fuel  deple- 
tion steering.  The  5 rad/sec  autopilot  has  adequate 
responsiveness  to  follow  this  probile,  as  well  as  the 
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desired  noise  attenuation  characteristics  at  higher  fre- 
quencies. The  optimized  S rad/sec  autipilot  has  approximately 
the  same  response  characteristics  as  the  8 rad/sec  system, 
but,  of  course,  has  a narrower  bandwidth.  This  is  an  impor- 
ant  advantage,  as  was  pointed  out  in  IV. 

Optimized  Design 

I.  The  optimized  5 rad/sec  design  gave  much  improved  performance 
over  the  standard  5 rad/sec  autopilot.  The  shape  of  the  fuel 
depletion  arc  for  this  configuration  is  very  much  like  that 
of  the  8 rad/sec  autopilot  in  that  the  switch  to  steering 
occurs  at  about  50  seconds,  rather  than  at  45  seconds. 

Delaying  the  switch  back  to  steering  to  about  TgQ  = 10 
seconds  (or  t = 50  seconds)  lowers  the  required  turning 
rate  during  the  fuel  depletion  maneuver  by  stretching  the 
maneuver  over  a longer  period.  This  also  seems  to  be  about 
the  optimal  time  to  switch  back  to  steering  in  order  to 
have  the  maximum  time  to  null  the  V errors. 

g 

The  numerical  comparison  of  the  terminal  conditions  for  the 

OFFSETT  and  STEER  simulation  runs  is  given  in  Table  21.  It  should  be 

noted  that  the  residual  Vg  values  shown  in  this  table  could  in  most 

cases  be  reduced  by  performing  a fine  countdown  during  the  last  few 

seconds  of  the  burn.  The  V and  V _ values  in  the  table  are  the 

g cap 

results  of  shutting  down  the  engine  at  the  next  autopilot  cycle  after 

either  1)  the  x-component  of  the  vector  goes  to  zero  (V^  steering) 

or  2)  the  magnitude  of  V goes  to  zero  (fuel  depletion  steering) . 

c»p 

Since  the  rate  of  change  of  V is  equal  to  the  thrust  acceleration 

cap  2 

and  for  the  constant  thrust  case,  the  acceleration  is  about  300  ft/sec 

at  the  end  of  the  burn,  we  can  see  that  Vcap  changes  by  about  10  fps 

between  autopilot  cycles.  Obviously,  further  refinements  can  be  made. 

The  lowest  residual  Vg  error  (1.24  fps)  was  achieved  by  the  optimized 

5 rad/sec  autopilot  (Run  1/23) . This  run  also  ended  with  the  lowest 

terminal  pitch  rate  (0.003  deg/sec). 
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Table  21.  Comparison  of  STEER  and  OFFSETT 
terminal  conditions. 


Run  # 

V 

g 

V 

gx 

V 

gz 

V 

cap 

§ 

Autopilot/ 
Steer ing/TFR 

Mod 

1/01 

5.13 

-4.37 

-2.68 

-4.9 

0.000 

5/01/0 

V9 

1/02 

5.11 

-4.35 

-2.67 

-4.9 

0.161 

5/01/0 

Vc 

1/03 

1.60 

-1.53 

-0.46 

2259.2 

-0.250 

5/01/0 

V 

c 

1/04 

1.60 

-1.53 

-0.46 

2259.2 

-0.225 

5/01/5 

V 

c 

1/05 

7.67 

-4.19 

-6.42 

-0.8 

2.092 

5/01/0 

FI 

1/06 

19.46 

-16.63 

-10.11 

-3.9 

6.797 

5/02/0 

FI 

1/07 

24.58 

-20.31 

-13.84 

-3.9 

1.222 

5/02/2 

FI 

1/08 

24.60 

-20.37 

-13.80 

-4.0 

1.166 

5/02/2 

FI 

1/09 

26.81 

-19.90 

-17.96 

0.2 

-1.130 

5/02/2 

FI 

1/10 

11.61 

-9.30 

-6.94 

-4.6 

0.307 

5/01/2 

FI 

1/11 

15.09 

-14.51 

-4.15 

-7.5 

0.239 

5/01/2 

FD 

1/12 

58.95 

-57.82 

-11.50 

-2.5 

-2.109 

5/01/2 

FD 

1/13 

397.94 

44.41 

-395.46 

365.8 

4.177 

2/01/2 

FI 

1/14 

8.18 

-6.06 

5.50 

-5.6 

-1.903 

2/02/2 

FI 

1/15 

2.22 

-1.84 

-1.24 

-1.0 

0.082 

8/01/2 

FI 

1/16 

19.77 

-17.00 

-10.08 

-3.8 

1.921 

8/02/2 

FI 

1/17 

2.85 

-2.83 

0.29 

-2.6 

0.039 

5*/01/2 

FI 

1/18 

11.88 

-11.22 

-3.90 

-4.2 

1.787 

5*/02/2 

FI 

1/19 

12.07 

-11.35 

-4.10 

-4.1 

1.781 

5*/02/2 

FI 

1/20 

1.04 

-0.69 

-0.77 

0.2 

-1.958 

5V02/2 

FI 

1/21 

72.12 

59.55 

-36.32 

-5.2 

0.141 

8/01/2 

FI 

1/22 

76.40 

63.37 

-38.47 

-4.5 

0.115 

8/01/2 

FI 

1/23 

1.24 

0.91 

-0.84 

0.0 

0.003 

5*/01/2 

F 

OFFSETT  RUNS 

2/100  0.01  0.01  0.005  0.0  N/A  (Not  applicable)  V( 

.26  0.1 

.75  0.2 


APPENDIX  A1 


I I 


VEHICLE  ROTATIONAL  TRANSFER  FUNCTION 


Rotational  Dynamics  (Pitch  or  Yaw  Axis) 


Using  Eulers'  Equation 


18  = DM  = TJt  6 + M l l 6 + 16  (See  Table  1) 

eg  e e eg  e 


assuming 


(sin  6=6) 


The  first  term  on  the  right-hand  side  is  the  torque  about  the 
vehicle  eg  caused  by  the  defection  of  the  thrust  line;  the  second  term 
is  the  inertial  reaction  torque  due  to  lateral  translation  of  the  engine 
nozzle;  and  the  third  term  is  due  to  the  angular  acceleration  of  the 
engine.  Taking  the  Laplace  transform  of  both  sides 


is“e(s)  = 


TAcg«(s>  + CVe*cg  + Ie)S  6(s) 


9 (S) 
6 (s) 


Tl 


-£2. 


“twd2 


+ 1 


with 


Til 


TWD 


_C£_ 


Mil  + L 
e e eg  e 


If  we  substitute  typical  vehicle  parameters  into  this  expression,  we 
find  that  = 80  rad/sec.  This  frequency  is  much  higher  than  the 

autopilot  crossover  frequency  and  for  the  purposes  of  this  study,  this 
effect  will  be  neglected. 
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APPENDIX  A2 


ACTUATOR  TRANSFER  FUNCTION 


The  actuator  open-loop  transfer  function  consists  of  A DC  gain 
factor,  a first-order  lag  with  time  constant,  t,  and  an  integrator, 
shown  in  Figure  A2-1.  With  the  unity  feedback  loop  added,  the  closed- 
loop  transfer  function  becomes 


Gfi<s) 


K„/t 


We  choose  t = 1/Kg  so  that  the  damping  ratio,  c is  equal  to  0.5  and  the 
natural  frequency,  un,  is  equal  to  Kg.  Then  the  closed  loop  transfer 
function  can  be  written  as 


G5(s) 


2 2 
s*  + K s + K 
s s 


I 


Figure  A2-1.  Actuator  Model. 


APPENDIX  A3 


DERIVATION  OF  X(s)/e(s) 


The  autopilot  and  steering  loop  models  described  in  Section  2.3 
employ  the  variable  X to  represent  the  angular  position  of  the  thrust 
acceleration  in  the  pitch  or  yaw  plane  relative  to  the  initial  body  x- 
axis  that  would  be  sensed  by  the  zero-quantization  accelerometers  of 
an  idealized  IMU.  If  the  vehicle  rotates  about  its  pitch  (y)  or  yaw 
(z)  axis  with  an  angular  velocity  § the  IMU-sensed  acceleration  will 
have  an  x-axis  component 


= T cos  6 /M  + l S' 
XIMU  e IMU 


and  a component  perpendicular  to  the  x-axis 


iMcr 


T sin  6e/M  + *IMU  0 


(A-l) 


(A-2) 


where  6£  is  the  angular  displacement  of  the  engine  from  the  x-axis 

(defined  positive  in  the  direction  that  produces  a positive  0)  and  iIMU 

is  the  distance  of  the  IMU  forward  of  the  vehicle  center  of  mass. 

We  neglect  the  centripetal  acceleration  term,  and  assume 

that  6g  is  small  enough  that 

cos  6e  = 1 (A-3) 

sin  6e  = 60  (A-4) 


The  body-axis  components  of  the  IMU-sensed  acceleration  are  then 
approximated  by 


aY  » T/M 
aimu 


IMU 


- t/m  <se  + *IMUe 


(A- 5) 

(A-6) 
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UTiltI  <,e  ax  ' tlie  angle  of  the  IMU-sensed  acceleration 


Assuming  that  a. 

'IMU  “imu 
relative  to  the  body  x-axis  is  given  by  the  ratio 


’IMU 


IMU 


a + e 

6e  + ^7ST  9 


(A-7) 


The  angle  X of  this  IMU-sensed  acceleration,  measured  relative  to 
the  same  inertially  fixed  coordinates  in  which  0 is  measured,  must  be 

given  by  the  sum  6 + a /a  , whence 

IMU  XIMU 


1 " 6 - SE  + 1IMO/T/M| 


[1IMO/T/M|  «] 

The  transfer  function  A(s)/8(s)  may  therefore  be  expressed  as 

~ [ i + s2 

' L 1 3 Tt7m ) . 


(A-8) 


* <«) 

0(s) 


6F(s) 

9(3) 


(A-9) 


But 


10  * 


T l 

e e 


(A -10) 


where 


therefore 


I = Moment  of  inertia  about  the  axis 

of  rotation  for  which  8 is  defined 

»e  = Moment  arm  between  the  engine 
hinge  point  and  the  vehicle 
center  of  mass 


and 


5e(s) 


- s2I 
0(S)  “ TJl 

e 


Ms) 
0 (a) 


1 + s' 


IMU 


(T/M) 


I 

11. 


(A-ll) 


(A-12) 
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APPENDIX  A4 

j 

RELATIONSHIP  BETWEEN  X,  XAy,  AND  Xy  ] 

g -j 


We  suppose  that  Figure  5 in  Section  2.3  represents  transfer 
functions  in  a pitch  plane  whose  inertially  fixed  x-axis  is  the 
reference  for  the  previously  defined  6 and  X.  If  we  further  assume 
small  engine  deflections  such  that  the  thrust  acceleration  vector  is 
approximately  parallel  to  the  vehicle  x-axis  we  have  from  Figure  A4-1. 

mTs 

AVz(mV  = I / Xdt 

(m-l)Ts 

and 


AVX (mTg) 


T 

M 


T3 


Z 


Figure  A4-1.  Definition  of  angles  in  inertial  frame. 


Therefore,  using  the  fact  that  tan  x = x for  small  x 

mT_ 


AV  (mT  ) 


£ / 

(m-1) T 
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The  integral  is  equivalent  to 


XAV(mTs> 


(m-l)T 

/ 


s 


Xdt 


o 


which  can  be  represented  pictorially  by  Figure  A4-2.' 


XU) 


Figure  A4-2.  Relation  between  X(s)  and  X^v(s). 


The  representation  of  the  effects  of  Vg  in  Figure  5 is  based  on  the 
assumption  that 


and  the  assumption  that  the  initial  direction  of  Vg  points  along  the 
x-axis  of  the  pitch  (or  yaw)  plane.  For  the  case  of  pitch  plane  motion, 
it  is  assumed  in  Figure  5 in  Section  2.2  that 


V 

x 


T 

M t 


V 

z 


t 
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and  also  that 


V 


gx 


go 


V 


gz 


■ -If  ldt 


Finally,  is  is  assumed  that 


whence 


/ 


Adt 


Thus,  the  effect  of  sampling  to  determine  Ay  every  Ts  seconds  can 
be  represented  by  the  following  diagram.  g 


Figure  A4-3.  Relationship  between  Ay  (s)  and  A(s). 

g 

In  Figures  5 and  6,  Av  (s)  has  been  replaced  by  the  more  general  term 

g 

AUSF(s)  and  the  term  -1/Tg0  has  been  changed  to  -Y/Tgo  since  the  Y/Tg 
term  represents  the  1/T  sensitivity  for  Vg  steering  and  4/Tgo  sen- 
sitivity of  fuel  depletion  steering  as  discussed  in  Section  3.3. 
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Figure  A5-1.  Geometry  of  the  fuel  depletion  arc. 


In  this  figure,  0 represents  the  fuel  depletion  angle,  the 
chord  length  c represents  the  magnitude  of  Vg  and  the  arc  length  S 
represents  Vcap. 

Note  that 
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APPENDIX  A6 


CALCULATION  OF  SPECIFIC 
FLIGHT  CONTROL  SYSTEM  PARAMETERS 


The  relationship  between  the  approximate  autopilot  crossover 
frequency  and  w'  is 


uxo  " “ “2 


to-  — to)  /to)  ’ 

2 xo 


Once  0)2  is  determined,  W2  is  found  from 


w2  = 2 “2 


Since  we  have  selected  = 0.10 


— 0.10  uij 


The  w-plane  pole  of  the  compensator  is  given  by 

T/2 

W3  T, 


and  T. 


($'/to)2)  + TDAp/2  + Tact. 


We  have  selected  * 0.25  and 

T =0.03  and 
DAP 


0.020833 


Finally  the  compensator  gain  is 


where 


and 


Using  these  equations  and  values  for  the  fixed  parameters,  we  generate 
Table  7.  With  the  substitution 


for  this  simple  compensator,  we  can  determine  the  z-plane  coefficients 
from  the  equations 


i 
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APPENDIX  B 

COORDINATE  TRANSFORMATIONS 


f 

i 

j 

i 
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The  coordinate  transformations  required  for  the  simulation 

programs  are  presented  herein.  The  relative  orientation  of  the  body 

frame  and  the  inertial  frame  is  shown  in  Figure  B-l,  and  it  assumes 

that  the  platform  axes  are  initially  aligned  with  the  inertial  axes. 

* ? 

The  platform-to-body  transformation  matrix,  PB,  and  its  inverse,  the 

* i 

body-to-platform  matrix,  BP,  are  used  to  transform  inertially-ref erred  j 

vectors  into  the  body  frame  and  body  referenced  vectors  into  the  inertial  j 

(platform)  frame.  These  transformation  matrices  are  updated  at  the  J 

autopilot  loop  cycle  time.  j 

BODY  FRAME  TO  INERTIAL  ] 

FRAME  TRANSFORMATION  3 


Zii 


Figure  B-l.  Body  frame  to  inertial  frame  transformation. 
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The  engine  frame-to-body  frame  transformation  matrix  is  presented 

in  Figure  B-2.  It  is  used  to  transform  the  engine  orientation  angles 

into  a body  referenced  vector,  ENG,  which  is  then  transformed  into  the 

inertial  frame  by  multiplying  it  by  the  body-to-platform  transformation 
* — 
matrix,  BP.  This  inertially  referred  vector,  ENGPLAT,  represents  the 

direction  of  the  engine  thrust  vector. 

THRUST  VECTOR 


BODY  LONGITUDINAL  AXIS 


ENGINE  CENTERLINE 


The  direction  of  the  thrust  vector  is  referenced  to 
the  vehicle  body  axes  using  the  ENG  vector  where 

cos  6y  cos  6 

P 

ENG  = -sin  6y  cos  6 

P 

sin  6_ 


Figure  B-2.  Engine  frame  to  body  frame  transformation. 


APPENDIX  C 


SIMULATION  DESCRIPTION 


To  evaluate  the  design  of  the  guidance,  steering,  control  and 
attitude  measurement  systems,  digital  computer  programs  were  developed 
and  run  to  simulate  these  systems  independently,  as  well  as  in  an 
integrated  system.  Three  major  programs  (OFFSETT,  STAR,  and  STEER) , 
and  five  subroutines  (PERT,  AIRSBAE,  AIRSME2,  AIRSME4 , and  AIRSINC) 
were  developed  or  used  in  this  study. 

OFFSETT  is  a simulation  of  the  guidance  system  with  simplified 
vehicle  dynamics  operating  in  orbit  around  a spherical,  nonrotating 
earth.  STAR  simulates  the  steering  and  control  loops  and  the  vehicle 
and  actuator  dynamics,  but  has  no  guidance  routine  and  operates  in 
inertial  space.  These  two  major  programs  are  merged  into  STEER,  which 
is  a complete  simulation  of  the  integrated  guidance,  steering,  auto- 
pilot and  actuator  systems  coupled  with  the  complete  vehicle  dynamics 
in  earth  orbital  flight.  PERT  is  a subroutine  to  OFFSETT  that  is  used 
to  compute  sensitivities  to  velocity  errors.  The  AIRS  programs  are 
variations  of  the  program  that  simulate  the  attitude  data  processing. 

A description  of  each  of  these  programs  is  presented  below. 


OFFSETT 

The  OFFSETT  program  simulates  the  guidance  system.  The  vehicle 
is  represented  by  a point  mass,  with  no  steering  or  autopilot  inter- 
actions or  rotational  dynamics.  Normally,  the  thrust  and  the  mass 
flow  rate  are  constant,  but  the  same  thrust  model  (see  Appendix  D) 
used  in  STEER  is  available  for  use. 

The  program  inputs  are  the  initial  and  final  radius  vectors,  the 
time  of  flight,  the  initial  velocity,  and  the  reentry  angle  control 
tolerance.  The  program  begins  by  computing  an  impulsive  velocity  to  be 
gained  by  solving  the  Lambert  problem  using  the  initial  and  final 
radius  vectors  and  the  time  of  flight.  Then,  at  each  guidance  cycle, 
a position  offset  is  calculated  to  account  for  the  nonimpulsive 


nature  of  the  burn.  The  new  velocity-to-be-gained  is  calculated  from 
this  offset  position.  The  vehicle  may  be  steered  in  either  Vg  steering 
or  fuel  depletion  steering.  If  fuel  depletion  steering  is  used,  the 
reentry  angle  may  be  controlled  by  modulating  Vcap  to  change  the  fuel 
depletion  angle  and  hence,  change  the  pitch  profile  to  control  the 
reentry  angle,  while  at  the  same  time  maintaining  a fixed  time  of  flight. 
The  vehicle  dynamics  are  simulated  in  a differential  equation  loop, 
and  all  measurements  are  considered  ideal. 

The  OPFSETT  program  and  its  variations  are  used  as  an  "ideal” 
case  to  show  just  the  effects  of  the  guidance  computation  without 
steering  or  vehicle  lags. 


0FFSETT1 

This  program  is  identical  to  OFFSETT,  except  that  the  thrust  model 
explained  in  Appendix  D is  included  to  simulate  the  effects  of  thrust 
tailoff . 


PERT  is  a subroutine  to  OFFSETT  (and  OFFSETTl ) that  computes  the 
change  in  fuel  depletion  angle  caused  by  a small  variation  in  the 
velocity- to-be-gained.  The  PERT  subroutine  is  called  after  the 
guidance  calculations  at  a particular  time  have  been  made.  A small 
variation  in  the  vehicle  velocity,  and  hence  the  velocity-to-be-gained, 
is  made  and  the  perturbed  parameters  are  calculated  at  the  same  time  as 
the  unperturbed  parameters.  After  these  parameters  are  calculated,  the 
unperturbed  parameters  are  restored  and  the  simulation  continues. 

Sample  output  for  the  OFFSETT  program  with  the  PERT  feature 
activated  is  shown  in  Table  C-l. 


The  STAR  program  simulates  the  steering  loops,  autopilot  filter, 
engine  actuator,  attitude  measurement  and  processing  system  (using  the 
AIRS  subroutines)  without  inputs  from  the  guidance  program.  It  was 
used  to  evaluate  the  response  of  the  vehicle  to  test  steering  inputs  to 
examine  the  steering  concepts,  autopilot  filter  parameters  and  two 
different  methods  of  processing  the  attitude  data.  In  the  STAR  program, 
the  vehicle  may  be  steered  using  vg  steering  or  a test  steering  input 
may  be  used.  For  the  purpose  of  the  AIRS  study  presented  in  Section  4, 
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a combination  of  a parabolically  increasing  constant,  and  linearly  de- 
creasing rate  input  was  used.  In  the  STAR  program,  the  mass,  moments 
of  inertia,  and  center  of  gravity  location  are  linearly  varying  func- 
tions of  time.  The  thrust  is  constant  but  the  autopilot  filter  gain 
is  changed  as  a function  of  thrust  acceleration  to  maintain  a constant 
value  of  open  loop  gain. 

Sample  output  of  the  STAR  program  is  given  in  Table  C-2. 

AIRS 

The  variations  of  the  AIRS  programs  differ  in  the  methods  that 
are  used  to  extrapolate  and  interpolate  the  AIRS  band  angles  and 
compute  the  transformation  matrices.  In  general,  though,  the  AIRS 
programs  receive  the  platform-to-body  transformation  matrix  and  the 
values  of  receiver  and  driver  band  noise  from  the  main  program  and 
compute  quantization  and  deterministic  (table  look-up)  errors  to 
determine  the  measured  incremental  attitude  vector  and  the  measured 
body-to-platform  transformation  matrix.  The  variations  of  the  AIRS 
program  are  discussed  further  in  Section  4. 

STEER 

The  STEER  program  combines  all  aspects  of  the  previously  described 
programs  into  one  complete  simulation  program,  as  shown  in  Figure  C-l, 
the  simplified  STEER  flowchart.  Referring  to  Figure  C-l,  after 
initialization  of  the  parameters,  the  program  computes  the  guidance 
parameters  under  GUID.  These  parameters  include,  the  predicted  engine 
cut-off  time,  which  is  either  the  time  when  the  velocity-to-be-gained 
should  be  zero  or  when  the  V is  exhausted,  the  position  offset,  the 

C9p 

fuel  depletion  angle  if  fuel  depletion  guidance  is  used,  the  current 
velocity- to-be-gained,  and  the  reentry  angle.  Within  this  loop,  a 
solution  convergence  check  recomputes  the  parameters  until  the  position 
offset  and  guidance  solution  converges.  After  the  guidance  parameters 
are  computed,  the  program  enters  the  STEER  routine  > here  the  velocity- 
to-be-gained,  the  steering  error  vector,  and  the  rt.^e  command  vector 
are  computed.  The  rate  vector  is  multiplied  by  TD;>r,  to  compute  the 
incremental  attitude  command  vector.  The  program  tnen  enters  CONTRL 
where  this  vector  is  compared  to  the  incremental  attitude  feedback  and 
an  attitude  error  vector  is  computed.  This  error  signal  is  used  by  the 
autopilot  filter  to  generate  actuator  deflection  commands.  Under 
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AIRS  ANGLE  INCREMENT  s 0.5761  0.6006  0.5958 

TRUE  ANGLE  INCREMENT  =-0.0000  0.6019  0.5935 


UPDATE,  the  vehicle  mass,  center  of  gravity,  and  moment  of  inertia  data 
is  updated  and  if  the  thrust  tailoff  option  is  selected,  the  thrust,  mass 
flow  rate,  and  exhaust  velocity  is  modified.  The  vehicle  acceleration 
is  then  calculated.  The  actuator  commands  previously  computed  are  then 
used  in  the  ACTUATOR  subroutine  which  is  the  start  of  the  differential 
equation  routine  to  compute  the  actuator  response.  Under  DYNAM,  the 
actuator  response  and  the  rest  of  the  vehicle  parameters  are  used  to 
compute  the  vehicle  response.  The  program  enters  the  AIRS  subroutine 
where  the  actual  and  AIRS  program  computed  attitude  data  is  calculated. 

If  the  elapsed  burn  time  is  less  than  the  cut  off  time,  the  program 
returns  to  GUID.  However,  since  there  are  nominally  45  CONTROL  cycles 
between  GUID  cycles,  the  program  jumps  to  STEER.  There  are  nominally 
j 15  CONTRL  cycles  between  STEER  cycles,  so  the  program  then  jumps  to 

' CONTRL  and  exercises  that  loop.  This  process  continues  for  15  CONTRL 

! cycles  and  then  the  STEER  routine  is  exercised.  After  a total  of  45 

CONTRL  cycles  (and  3 STEER  cycles) , the  guidance  parameters  are  re- 
computed. The  program  continues  in  this  fashion  until  T = Tco  at  which 
j time  the  program  ends.  Sample  output  for  the  STEER  program  is  shown 

\ in  Table  C-3. 
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APPENDIX  D 


THRUST  TAILOFF  MODEL 


For  the  purposes  of  much  of  this  study,  the  rocket  thrust  was 
assumed  to  be  constant.  The  mass  flow  rate  was  also  constant,  and 
hence  the  mass  decreased  linearly  with  burn  time.  The  magnitude  of 
the  thrust  acceleration  was  given  by 


. _ Thrust 

T M - mt 
o 


The  AV  capability  of  the  vehicle  is  given  by  the  integral  of  the 
thrust  acceleration 

t t 


(Dl) 


/v*  = / 


cap  •/  Ardt 

F o o M - mt 


dt 


(D2) 


With  the  assumptions  stated  above,  this  integral  can  be  evaluated 
analytically  as 

Vcap  = T ln(M0/Mburn0ut)  (D3) 

m 

T 

Recognizing  — as  the  exhaust  velocity,  V , the  expression  becomes 
m 


cap 


vexln  ^o^burnout^ 


(D4) 


While  the  assumption  of  constant  thrust  and  linear  mass  change  is  a 
useful  simplification  in  the  design  process,  most  real  solid-fuel 
rocket  engines  exhibit  thrust  tailoff  characteristics,  marked  by 
decreasing  thrust,  mass  flow  rate,  and  exhaust  velocity  in  the  latter 
part  of  the  burn,  as  well  as  thrust  variations  throughout  the  burn, 
caused  by  unequal  burning  of  the  fuel. 

As  has  been  demonstrated  in  Section  2,  the  steering/autopilot/ 
vehicle  combination  is  stable  only  for  a range  of  value  of  open-loop  gains. 
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Since  the  gain  contribution  of  the  vehicle  is  TACg/I,  as  the  thrust 
decreases  rapidly  at  the  end  of  the  burn,  the  By stem  tends  toward  in- 
stability. This  effect  may  be  compensated  for  by  varying  the  steering 
and/or  autopilot  compensator  gain  as  a function  of  sensed  acceleration. 

To  correctly  represent  the  effects  of  thrust  tailoff,  it  is  not 
adequate  to  simply  decrease  the  thrust  to  zero  in  some  fashion.  The 
other  engine  and  vehicle  parameters  (m,  v„,  V___,  M,  Am)  must  be 

6X  Cap  J- 

simultaneously  changed  to  make  a consistent  set.  Since  the  AV  of 

cap 

the  engine  is  assumed  to  be  accurately  known  beforehand,  the  first 

requirement  of  consistency  is  to  vary  the  acceleration  such  that  the 

area  under  the  curve  (the  vehicle  AV  capability)  remains  the  same. 

Since  the  acceleration  decreases  with  the  thrust,  in  order  to  keep  the 

same  AV , the  burn  time  must  be  changed.  This  forces  a change  in  the 
cap 

mass  flow  rate  and  hence  the  exhaust  velocity  changes.  If  the  thrust 
and  mass  flow  rate  assumed  constant  up  to  a certain  time  into  the  burn 
the  integral  expression  for  AV  may  be  expressed  as 
±bo  tl  tbo 


AT(t)dt 


Mq  - mt 


tl  Mtl  - m(t)t 


If  we  take  T(t)  - T - K.  x and 

JL 


m(t)  ’ “max  " V 


where  t * (t  - tl) f from  tl  to  tbo. 


M (t)  = - (m  - K2t  )t 


Then  and  can  be  computed  from 


“bo  - Mtl  + mT 
3 


Where  x is  selected  to  make  Vcap  equal  the  desired  value.  The 
variation  of  thrust,  acceleration,  mass,  and  mass  flow  rate  are 
shown  in  Figures  Dl  through  D4,  Engine  Tailoff  Model. 
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gure  D-2.  Thrust 


APPENDIX  El 


EFFECT  OF  IMU  POSITION 
ON  ACCELERATION  AND  RATE  VECTORS  SENSED 

The  Inertial  Measurement  Unit  (IMU) , which  contains  the  acceler- 
ometer and  rate  gyros,  is  normally  located  in  the  post-boost  vehicle  or 
payload  vehicle  since  this  information  is  required  after  the  powered 
flight  portion  of  the  trajectory  is  over.  During  the  third  stage  burn, 
this  results  in  a large  moment  arm  between  the  center  of  gravity  of 
the  vehicle  and  the  IMU  location.  The  effect  of  this  moment  arm  is  to 
cause  errors  in  the  sensed  rate  and  acceleration  vectors  when  the 
vehicle  is  rotating.  These  errors  are  due  to  the  centripetal  and 
tangential  components  of  acceleration  caused  by  the  rotation.  Figure 
El  shows  this  pictorially. 


The  maximum  tangential  acceleration!  A^,^,  sensed  by  the  XMU  is 

the  maximum  angular  acceleration!  9 times  the  moment  armf  r. 

max 


Assuming  9 


20  deg/sec  , and  r - 10  ft.,  this  gives  a maximum 


tangential  acceleration  of  ATAN  * 3.495  ft/sec  • The  worst  case 

condition  occurs  when  the  magnitude  of  the  thrust  acceleration  A^, 

is  at  a minimum.  The  minimum  (ignition)  value  of  the  thrust  acceleration 
2 

is  A^  » 94.3  ft/sec  . Adding  these  twc  accelerations  vectorally  gives 
a worst  case  false  thrust  direction  indication  of  2.12  degrees.  Since 
this  is  a small  value,  and  since  the  maximum  value  of  9 occurs  early 
in  the  burn  (at  the  initiation  of  fuel  depletion  steering)  and  will  be 
corrected  by  guidance,  this  effect  will  be  neglected  for  the  purposes 
of  this  study. 

The  maximum  centripetal  acceleration,  A , sensed  by  the  IMU 
is  the  square  of  the  maximum  angular  rate,  9 , times  the  moment  arm,  r. 


acent  “ r 


Assuming  a maximum  angular  rotation  rate,  § = 25  deg/sec,  and 

again  r * 10  ft.,  this  gives  a maximum  centripetal  acceleration  of 
2 

1.904  ft/sec  . Since  this  value  is  small  compared  to  the  magnitude 
of  the  thrust  acceleration,  it  will  be  neglected  for  the  purposes  of 
this  study. 
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